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Abstract 



This monograph is derived from a series of six lectures I gave at the Centre de Recherches 
Mathematiques in Montreal, in March and June 2000, while titulary of the Aisenstadt Chair. 



Contents 



Precis 



1 Introduction 



2 Classical Equations 

2.1 Equations of motion 

2.2 A word on canonical formulations 

2.3 The irrotational case 

2.4 Nonvanishing vorticity and the Clebsch parameterization 

2.5 Some further remarks on the Clebsch parameterization . 



3 Specific Models 

[Tj Galileo -invariant nonrelativistic modefl . . . 

3.2 Lorentz-invariant relativistic model 

3.3 Some remarks on relativistic fluid mechanics 



4 Common Ancestry: The Nambu-Goto Action 

4.1 Light-cone parameterization 

4.2 Cartesian parameterization 

4.3 Hodographic transformation 

4.4 Interrelations- 



5 Supersymmetric Generalization 

|5.1 Chaply" gin gas with Grassmann variables] 

5.2 Super symmetry 

5.3 Supermembrane Connection 

5.4 Hodographic transformation 

5.5 Light-cone parameterization 

5.6 Further consequences of the supermembrane connection 



6 One-dimensional Case 

5.1 Specific solutions for the Chaplygin gas on a line 

|6.2 Aside on the integrability of the cubic potential in one dimension 

3.3 General solution for the Chaplygin gas on a line 

5.4 Born-Infeld model on a line] 

5.5 General solution of the Nambu-Goto theory for a (d=l)-brane (string) in two| 
spatial dimensions (on a plane) 



7 Towards a Non-Abelian Fluid Mechanics 

7.1 Proposal for non-Abelian fluid mechanics 

7.2 Non-Abelian Clebsch parameterization] . . 

7.3 Proposal for non-Abelian magnetohydrodynamicsj 



Solutions to Problems 



References 



Lectures on (Supersymmetric, Non-Abclian) Fluid Mechanics (and d-Branes) 



•5 



Precis 

During the March 2000 meeting of the Workshop on Strings, Duality, and Geometry in 
Montreal, Canada, I delivered three lectures on topics in fluid mechanics, while titulary of 
the Aisenstadt Chair. Three more lectures were presented in June 2000, during the Montreal 
Workshop on Integrable Models in Condensed Matter and Non-Equilibrium Physics. Here are 
brief descriptive remarks on the content of the lectures. 

1. Introduction - The motivation for the research is explained. 

2. Classical Equations - The classical theory is reviewed, but in a manner different from 
textbook discussions. 

(a) Equations of motion - Summary of conservation and Euler equations. 

(b) A word on canonical formulations - An advertisement of the method for finding the 
canonical structure for the above (developed with L.D. Faddeev). 

(c) The irrotational case - C. Eckart's Lagrangian and a relativistic generalization for 
vortex-free motion. 

(d) Nonvanishing vorticity and the Clebsch parameterization - In the presence of vor- 
ticity, the velocity Chern-Simons term (kinetic helicity) provides an obstruction to 
the construction of a Lagrangian for the motion. C.C. Lin's method overcomes the 
obstruction, and leads to the Clebsch parameterization for the velocity vector. 

(e) Some further remarks on the Clebsch parameterization - Properties and peculiarities 
of this presentation for a 3-vector. 

3. Specific Models - Nonrelativistic and relativistic fluid mechanics in spatial dimensions 
greater than one. 

(a) Galileo-invariant nonrelativistic model - The Chaplygin gas [negative pressure, in- 
versely proportional to density] is studied, selected solutions are presented, unex- 
pected symmetries are identified. 

(b) Lorentz-invariant relativistic model - The scalar Born-Infeld model is found to be 
the relativistic generalization of the Chaplygin gas, and shares with it unexpected 
symmetries. 

(c) Some remarks on relativistic fluid mechanics - Dynamics for isentropic relativistic 
fluids is given a Lagrangian formulation, and the Born-Infeld model is fitted into 
that framework. 

4. Common Ancestry: The Nambu-Goto Action - Both the Chaplygin gas and the Born- 
Infeld model devolve from the parameterization-invariant Nambu-Goto action, when spe- 
cific parameterization is made. 

(a) Light-cone parameterization - Chaplygin gas is derived. 

(b) Cartesian parameterization - Born-Infeld model is derived. 
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(c) Hodographic transformation - Chaplygin gas is derived (again). 

(d) Interrelations - The Chaplygin gas and Born-Infeld are related because (1) the former 
is the nonrelativistic limit of the latter; (2) both descend from the same Nambu-Goto 
action. 

5. Supersymmetric Generalization - Fluid mechanics enhanced by supersymmetry. 

(a) Chaplygin gas with Grassmann variables - Vorticity is parameterized by Grassmann 
variables, which act like Gaussian potentials of the Clebsch parameterization. 

(b) Supersymmetry Supercharges, transformations generated by 
them, and their algebra. 

(c) Supermembrane connection - Supermembrane Lagrangian in three spatial dimen- 
sions. 

(d) Hodographic transformation - Supersymmetric Chaplygin gas in two spatial dimen- 
sions is derived. 

(e) Light-cone parameterization - Supersymmetric Chaplygin gas in two spatial dimen- 
sions is derived (again). 

(f ) Further consequences of the supermembrane connection - Hidden symmetries of the 
supersymmetric model. 

6. One-dimensional Case - The previous models in one spatial dimension are completely 
integrable. 

(a) Solutions for the Chaplygin gas on a line - Some special solutions are presented; 
infinite number of constants of motion is identified; Riemann coordinates are intro- 
duced and the fluid equations as well as constants of motion are expressed in terms 
of them. 

(b) Aside on the integrability of the cubic potential in one dimension - The one-dimen- 
sional problem with pressure oc (density) 3 possesses the SO(2,l) "Schrodinger sym- 
metry" and the equations of motion, in Riemann form, become free. 

(c) General solution of the Chaplygin gas on a line - Solution obtained by linearization. 

(d) Born-Infeld model on a line - When formulated in terms of its Riemann coordinates, 
it becomes trivially equivalent to the Chaplygin gas. 

(e) General solution of the Nambu-Goto theory for a (d = l)-brane (string) in two 
spatial dimensions (on a plane) - The explicit string solution is transformed by a 
hodographic transformation to the Chaplygin gas solution, and a relation is estab- 
lished between this solution and the one found by linearization. 

7. Towards a Non-Abelian Fluid Mechanics - Motivation for this theory is given. 

(a) Proposal for non-Abelian fluid mechanics - A Lagrangian is proposed; it involves a 
non-Abelian auxiliary field whose Chern-Simons density should be a total derivative. 
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(b) Non-Abelian Clebsch parameterization (or, casting the non-Abelian Chern-Simons 
density into total derivative form) - Total derivative form for the non-Abelian Chern- 
Simons density is found, thereby generalizing the Abelian Clebsch parameterization, 
which achieves a total derivative form for the Abelian density. 

(c) Proposal for non-Abelian magnetohydrodynamics - Our proposal, which generalizes 
the one in Section 7.1 to include a dynamical non-Abelian gauge field, reduces in 
the Abelian limit to conventional magnetohydrodynamics. 



s 
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1 Introduction 

Field theory, as developed by particle physicists in the last quarter century, has enjoyed a 
tremendous expansion in concepts and calculational possibilities. 

We learned about higher and unexpected symmetries, and discovered evidence for partial 
or complete integrability facilitated by these symmetries. We appreciated the relevance of 
topological ideas and structures, like solitons and instantons, and introduced new dynamical 
quantities, like the Chern-Simons terms in odd-dimensional gauge theories. We enlarged and 
unified numerous degrees of freedom by introducing organizing principles such as non-Abelian 
symmetries and super symmetries. Indeed, application of field theory to particle physics has 
now been replaced by the study of fundamentally extended structures like strings and mem- 
branes, which bring with them new mathematically intricate ideas. 

Thinking about research possibilities, I decided to investigate whether the novelties that 
we have introduced into particle physics field theory can be used in a different, non-particle 
physics, yet still field-theoretic context. In these Aisenstadt lectures I shall describe an ap- 
proach to fluid mechanics, which is an ancient field theory, but which can be enhanced by the 
ideas that we gleaned from particle physics. 

As an introduction, I shall begin with a review of the classical theory. Mostly, I duplicate 
what can be found in textbooks, but perhaps the emphasis will be new and different. After 
this I shall describe how some instances of the classical theory are related to d-branes and 
how this relation explains some integrability properties of various models. I shall then show 
how the degrees of freedom can be enlarged to accomodate supersymmetry and non-Abelian 
structures in fluid mechanics. A few problems are scattered throughout; solutions are given at 
the end of the text, before the references. 

New work that I shall describe here was done in collaboration with D. Bazeia, V.P. Nair, 
S.-Y. Pi, and A. P. Polychronakos. Textbooks for the classical theory, which I recommend, are 
by Landau and Lifschitz ||] as well as by Arnold and Khesin p. 
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2 Classical Equations 
2.1 Equations of motion 

We begin with nonrelativistic equations that govern a matter density field p(t, r) and a velocity 
field vector v(t,r), taken in any number of dimensions. The equations of motion comprise a 
continuity equation, 

J^(t,r) + V-(p(t,r)i>(t,r))=0 (1) 

which ensures matter conservation, that is, time independence, of N = J dr p, and Euler's 
equation, which is the expression of a nonrelativistic force law. 

r) + v(t, r) • Vw(t, r) = f(t, r) (2) 

Here pv is the current j and / is the force. We shall deal with an isentropic fluid, that is, 
entropy is constant and does not appear in our theory. Also we ignore dissipation and take 
the force to be given by the pressure P: f = — i VP. For isentropic motion P is a function 
only of p, so / can also be written as —'VV'(p): 

f = —-VP = —W'(p) (3) 
P 

with the dash (also known as "prime") designating the derivative with respect to argument. 
V'(p) is the enthalpy, pV'(p) — V(p) = P(p), and y 7 P' (p) = s is the speed of sound. (Those 
familiar with the subject will recognize that I am using an Eulerian rather than a Lagrangian 
description of a fluid j|.) 

The dynamics summarized in (|l|) and (0) and the definition (S) may be presented as 
continuity equations for an energy momentum tensor. The energy density £ = T°° 



together with the energy flux 



obey 



£ = \pv 2 + V(p) = T°° (4a) 



T i° = pv i{\v 2 + V) (4b) 



JLr°° + djT j0 = 0. (4c) 



Similarly the momentum density, which in the nonrelativistic theory coincides with the current, 



V { = pv { = T oi (5a) 
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and the stress tensor T %3 

T ij = S ij ( P V' -V)+ pv i v j = 5 ij P + pv i v j (5b) 

satisfy 

+ djT'' = 0. (5c) 

at 

Note that T m ^ T to because the theory is not Lorentz invariant, but T lJ = T Jt because it is 
invariant against spatial rotations. [Thus T^ v is not, properly speaking, a "tensor", but an 
energy-momentum "complex" .] 

A simplification occurs for the irrotational case when the vorticity 

Wij = d{V° — djV 1 (6) 

vanishes. For then the velocity can be given in terms of a velocity potential 9, 

v = V6 (7) 

and equation (|2|) can be replaced by Bernoulli's equation. 

The gradient of (§) gives ©, with help of (|) and (@). 



Problem 1 In the free Schrodinger equation for a unit-mass particle, ifr^% = —\V 2 ijj, set 
ip = p l l 2 e %e l h , and separate real and imaginary parts. Show that the resulting equations are 
like those of fluid mechanics. What is the velocity? Is vorticity supported? What is the 
force /? 



Our equations can be presented in any dimensionality, but we shall mostly consider the 
cases of three, two, and one spatial dimensions. In the first case, the vorticity is a (pseudo-) 
vector 

u) = Vxu (9) 

in the second, it is a (pseudo-)scalar 

uj = £ ij diV j (10) 

while the last, lineal case is always simple because there is no vorticity and the velocity can 
always be written as the derivative (with respect to the single spatial variable) of a potential. 
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Dynamics of any particular system is most economically presented when a canonical/ action 
formulation is available. To this end we note that the above equations of motion can be 
obtained by (Poisson) bracketing with the Hamiltonian 



H = j dr (\pv 2 + V(p)) = J dr£ 



provided the nonvanishing brackets of the fundamental (p, v) variables are taken to be 

{v i (r), P (r')} = d i S(r-r') 
u)ij(r) 



(11) 
(12) 
(13) 



{v l {r),v J {r')} 



p{r) 



8(r 



(14) 



(The fields in curly brackets are at equal times, hence the time argument is suppressed.) An 
equivalent, more transparent version of the algebra (14) is satisfied by the field momentum 
density, 



V = pv . 

As a consequence of (|l4|) we have 

{V i (r),p(r')} = p(r)d i 5(r-r') 
{P l {r),V j {r')} = V j {r)di5{r - r') +V l (r')d j 5(r - r'). 



(15) 



(16) 



This is the familiar algebra of momentum densities. One verifies that the Jacobi identity is 
satisfied @. 

Naturally one asks whether there exists a Lagrangian whose canonical variables lead to the 
Poisson brackets (|i~4|) or (16) and to the Hamiltonian (11). In more mathematical language, 
we seek a canonical 1-form and a symplectic 2-form that lead to the algebra (O) or (16). 



Problem 2 Second-quantized Schrodinger fields satisfy equal-time commutation (anticom- 
mutation) relations, when describing bosons (fermions): \tp(r), ip*(r')\± = S(r — r'). Show 
that the algebra ( |l~6|) is reproduced (apart from factors of ifi) when p = 'P = Im.h^'Vtp. 
Since in the nonrelativistic theory "P = j, find j in terms of p and v, with v as determined in 
Problem |]. 



2.2 A word on canonical formulations 

I shall now describe an approach to canonical formulations of dynamics, publicized by Faddeev 
and me ||], which circumvents and simplifies the more elaborate approach of Dirac. 
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We begin with a Lagrangian that is first order in time. This entails no loss of generality 
because all second-order Lagrangians can be converted to first order by the familiar Legendre 
transformation, which produces a Hamiltonian: H(p,q) = pq — L(q,q), where p = dL/dq 
(the over-dot designates the time derivative). The equations of motion gotten by taking 
the Euler-Lagrange derivative with respect to p and q of the Lagrangian L(p,p;q,q) = pq — 
H(p, q) coincide with the "usual" equations of motion obtained by taking the q Euler-Lagrange 
derivative of L(q,q). [In fact L(p,p;q,q) does not depend on p.] Moreover, some Lagrangians 
possess only a first-order formulation (for example, Lagrangians for the Schrodinger or Dirac 
fields; also the Klein-Gordon Lagrangian in light-cone coordinates is first order in the light-cone 
"time" derivative). 

Denoting all variables by the generic symbol £ l , the most general first-order Lagrangian is 

L = ai (OC-H(0. (17) 

Note that although we shall ultimately be interested in fields defined on space-time, for present 
didactic purposes it suffices to consider variables £ J (t) that are functions only of time. The 
Euler-Lagrange equation that is implied by ( |i~7| ) reads 



where 



M()e = (is) 

MO = ^ - ^ <*> 

The first term in (|17]) determines the canonical 1-form: Oj(£)£*di = aj(£)d£ J , while gives 
the symplectic 2-form: da;(£)df = §/y(£) df d£?. 

To set up a canonical formalism, we proceed directly. We do not make the frequently 
heard statement that "the canonical momenta dL/d^ 1 = Oj(^) are constrained to depend on 
the coordinates and we do not embark on Dirac's method for constrained systems. 

In fact, if the matrix fij possesses the inverse there are no constraints. Then (|T§|) 
implies 

C = f ij (0^-- (20) 

When one wants to express this equation of motion by bracketing with the Hamiltonian 
one is directly led to postulating the fundamental bracket as 

= -/%)• (22) 

The Poisson bracket between functions of £ is then defined by 

Wfl,«(0)-^^. (23) 
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One verifies that (^) satisfies Jacobi identity by virtue of (|l9|). 

When fij is singular and has no inverse, constraints do arise, and the development becomes 
more complicated (see ||). 

Our problem in connection with (jlj) and ( |l6|) is in fact the inverse of what I have here 
summarized. From (|14|) and (16), we know the form of f %J and that the Jacobi identity 
holds. We then wish to determine the inverse /$,-, and also aj from (|l9|). Since we know the 



Hamiltonian from (11), construction of the Lagrangian ( |17D should follow immediately. 

However, an obstacle may arise: If there exists a quantity C(£) whose Poisson bracket with 
all the vanishes, then 

o = {e,c(m = -r 3 ^jC(o . (24) 

That is, f l i has the zero mode ^jC(C)> an d the inverse to namely the symplectic 2-form 
fij, does not exist. In that case, something more has to be done, and we shall come back to 
this problem. 

Totally commuting quantities like C(£) are called "Casimir invariants". Since they Poisson- 
commute with all the dynamical variables, they commute with the Hamiltonian, and are 
constants of motion. But these constants do not reflect any symmetry of the specific Hamilto- 
nian, nor do they generate any infinitesimal transformation on £ l , since the {C(£),£*} bracket 
vanishes. 

As will be demonstrated below, the algebra (0), (16) admits Casimir invariants, which 
create an obstruction to the construction of a canonical formalism for fluid mechanics; this 
obstruction must be overcome to make progress. (In the Lagrangian formulation of fluid 
mechanics these Casimirs are related to a parameterization-invariance of that formalism 0].) 



2.3 The irrotational case 

We now return to the specific issue of determining the fluid dynamical Lagrangian. The 
problem of constructing a Lagrangian which leads to ([b3|) and ( [Tq) can be solved by inspection 
for the irrotational case, with vanishing vorticity [see (J6h]. For then the velocity commutator 
in fll4| ) vanishes and (|7j) shows that the first equation in (|14|) can be satisfied by taking p and 
9 to be canonically conjugate. 

{0(r),p(r')} = S(r-r') (25) 

Thus the Lagrangian reads 

^irrotational = J dr(9p - H v= ^ d ) (26) 

where H is given by (|ll]) with v taken as in (^). The form of this Lagrangian can be understood 
by the following argument, due to C. Eckart H. 
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Consider the Lagragian for ./V point-particles in free nonrelativistic motion. With the mass 
m set to unity, the Galileo- invariant, free Lagrangian is just the kinetic energy. 



N 



n=l 



In a continuum description, the particle-counting index n becomes the continuous variable r, 
and the particles are distributed with density p, so that Yl n =i v n(t) becomes J dr p(t, r)v 2 (t, r). 
But we also need to link the density with the current j = pv, so that the continuity equation 
holds. This can be enforced with the help of a Lagrange multiplier 6. We thus arrive at the 
free, continuum Lagrangian. 

; Galileo _ / j / 1 2 



Z Gahlco = J dr (^l pv 2 + + V • (pv))) (28) 

Since ^Galileo j g £ rg j. orc j er j n time and the canonical 1-form j drOp does not contain v, the 
latter may be varied, evaluated, and eliminated B. Doing this, we find 

pv - pV6 = (29) 

and we conclude that V# is the velocity or, more precisely, that V# is the v derivative of the 
kinetic energy, that is, the momentum p, which in this nonrelativistic setting coincides with v. 



V6 = ^-W =p = v (30) 

Substituting this in (|28|), we obtain 

L Galilco = I dr (e p _l_ p{ ye?) (31) 



which reproduces (|2q ) with the interaction V(p) in ( [lT|) set to zero, and leads to the free 
version of the Bernoulli equation of motion (^) . 

Taking the gradient gives 

v + v ■ Vv = 0. (33) 



Problem 3 The Lagrange density for the unit-mass Schrodinger equation can be taken as 
^Schrodinger = *^'0*^V' — ' What form does this take after V is represented by 

p\/2 e id/h'} Compare with (p6|). 
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Remarkably, the same equation ( |33| ) emerges for a kinetic energy T(v) that is an arbitrary 
function of v. This will be useful for us when we study a relativistic generalization of the 
theory. If we replace (28) by 

L = J dr +0(p + V-Gcw))) (34) 

and vary v to eliminate it, we get in a generalization of (|30|) 



dT 

W = ^p. (35) 

So in the general case, it is the momentum - the v derivative of T(v) - that is irrotational. 
The Lagrange density becomes 



L = J dr(0p - ph(p) p= v e ) (36) 
where h(p) is the Legendre transform of T{v). 

dh 

Hp) =vp- T(v) — = v (37) 



Again varying 6 in ( |36|) gives the continuity equation 

n SL ° ■ I' A 9h (p) fa 
U = ——- = p — / dr p — • — 



59 r J r dp 

p — I dr pv • — V( 



Varying p gives 



Sp 

Taking the gradient, this implies with the help of (|3 



p + V-(pv). (38) 



= ^ = -9 - h{p). (39) 



diQ = -v ■ —p 

or 1 

= -v J tt-P 
ori 

i dp 1 d , 



dv k dri V ^ ^ 



On the other hand, (|35| ) implies that 



dv k 



(41) 



The two are consistent, provided the free Euler equation holds, that is, 

v k + v j d jV k = (42) 
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(as long as dp 1 /dv k = d 2 T /dv l dv k has an inverse). 

Let me observe that free motion is here governed by a Lagrangian that is not quadratic 



and the free equations are not linear. Nevertheless, the equations of motion ( |38|) and (42) can 
be solved in terms of initial data. 

p(t = 0,r) = p (r) (43) 
v(t = 0,r) = v (r) (44) 

Upon determining the retarded position q(t, r) from the equation 

q + tv (q)=r (45) 

one verifies that the solution to the free equations reads 

v(t,r)=v {q) (46) 
i da 1 i 

p(t,r)=p Q (q)\det-^-\. (47) 

A final remark: Note that the free Bernoulli equation (||) coincides with the free Hamilton- 
Jacobi equation for the action. 

2.4 Nonvanishing vorticity and the Clebsch parameterization 

We now return to our original Galileo-invariant problem and enquire about the Lagrangian 
for velocity fields that are not irrotational, that is, whose vorticity is nonvanishing. Here we 
specify the spatial dimensionality to be 3, and observe that the algebra (14) possesses a zero 
mode, since the quantity 

C(v) = J dre^tfdjV* = J dr v ■ u (48) 

(Poisson) commutes with both p and v. So the symplectic 2-form does not exist: in the 
language developed above, f iJ has no inverse. (Notice that C, also called the "fluid helic- 
ity", coincides with the Abelian Chern-Simons term for v [Jj]].) (In the irrotational case with 
vanishing u, the obstacle obviously is absent.) 

To make progress, one must neuteralize the obstruction. This is achieved in the following 
manner, as was shown by C.C. Lin @. 

We use the Clebsch parameterization for the vector field v. Any three-dimensional vector, 
which involves three functions, can be presented as 

v = V9 + aV/3 (49) 

with three suitably chosen scalar functions 9, a, and (3. This is called the "Clebsch parameteri- 
zation", and (a, (5) are called "Gaussian potentials" ||. In this parameterization, the vorticity 
reads 



u: = Va x V/3 



(50) 
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and the Lagrangian is taken as 

L = - J drp(9 + a/3) - H v =vo+ a vp (51) 



with v in H expressed as in (f49|). Thus (p, 9) remain canonically conjugate but another 
canonical pair appears: (pa, (3). The phase space is 4-dimensional, corresponding to the four 
observables p and v, and a straightforward calculation shows that the Poisson brackets (|14|) 
are reproduced, with v constructed by (fiS|). 

But how has the obstacle been overcome? Let us observe that in the Clebsch parameteri- 
zation C is given by 

C = j dre ijk di6d ja d k (3 (52) 

which is just a surface integral 

C = [ dS-(9u). (53) 



In this form, C has no bulk contribution, and presents no obstacle to constructing a symplectic 
2-form and a canonical 1-form in terms of 9, a and (3, which are defined in the bulk, that is, 
for all finite r. 

Lin gave an Eckart-type derivation of Q5l|): Return to ^Galileo j n (pgj^ anc j adc [ a father 
constraint, beyond the one enforcing current conservation 

^Galileo = j dr (^l pv 2 +d ^ + V . ( pv ^ _ pa + v . V/3) ^j (54) 



Setting the variation (with respect to v) to zero evaluates v as in p9|); eliminating v from (|54|) 
gives rise to (^l|). 

This procedure works in any number of dimensions, producing the same canonical 1-form 
in any dimension. This means that in two spatial dimensions, on the plane, where the (p, v) 
space is three-dimensional, the four-dimensional phase space (p, 9; pa, (5) is larger. Moreover, 
the analog to C in two spatial dimensions, that is, the obstruction to constructing a symplectic 
2-form, is not a single quantity: an infinite number of objects (Poisson) commute with p and v. 
These are 



Jdrpf^Y n = 0,±l,±2,.... (55) 



(Of course, the C n vanish in the irrotational case where there is no obstruction.) One can 
understand why there is an infinite number of obstructions by observing that phase space 
must be even dimensional, but (p, v) comprise three quantities on the plane. So a nonsingular 
symplectic form can be constructed either by increasing the number of canonical variables to 
four, or decreasing to two. The Lin/Clebsch method increases the variables. On the other 
hand, decreasing to two entails suppressing of one continuous and local degree of freedom, and 
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evidently this is equivalent to neutralizing the infinite number of global obstructions, namely, 
the C n . But I do not know how to effect such a suppression, so I remain with (|5l|). 

Note that ^Galileo - n (g) ; apart from a total time derivative, can also be written in any 
number of dimensions as 

^Galileo = J dr ( pT ( v ) _ p + a ft _ pv . ( V + aV/3)) 

= J dr( P T(v) - fid^e + adM (56) 

where we have introduced the current four-vector 

f = (cp,pv), (57) 

employed the four-vector gradient = (^Jj,V), and denoted the kinetic energy by T(v). 
These expressions will form our starting point for a relativistic generalization of the theory as 
well as a non-Abelian generalization. (That is why we have introduced the velocity of light in 
the above definitions; of course c disappears in the Galilean theory, as it has no role there.) 

Finally we observed that in one spatial dimension, where v can always be written as 9' 
and the vorticity vanishes, the irrotational canonical 1-form J dx9p is generally applicable 
and can equivalently be written as — 5 / dx dy p{x)e{x — y)v(y), where e is a ±l-step function, 
determined by the sign of its argument. Evidently this leads to a spatially nonlocal, but 
otherwise completely satisfactory canonical formulation of fluid motion on a line. 



2.5 Some further remarks on the Clebsch parameterization 

Let me elaborate on the Clebsch parameterization for a vector field, which was presented for 
the velocity vector in (f49|). Here I shall use the notation of electromagnetism and discuss 
the Clebsch parameterization of a vector potential A, which also leads to the magnetic field 
B = V x A. (Of course the same observations apply when the vector in question is the velocity 
field v, with V x v giving the vorticity.) 

The familiar parameterization of a three-component vector employs a scalar function 8 (the 
"gauge" or "longitudinal" part) and a two-component transverse vector At- A = V# + At, 
V • At = 0. This decomposition is unique and invertible (on a space with simple topology). 
In contrast, the Clebsch parameterization uses three scalar functions, 9, a, and /3, 

A = V9 + aV(3 (58) 
which are not uniquely determined by A (see below). The associated magnetic field reads 

B = Vax V/3. (59) 
Repeating the above in form notation, the 1-form A = A4 dr l is presented as 



A = d9+adf5 



(60) 
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and the 2-form is 

d^ = dad^. (61) 

Darboux's theorem ensures that the Clebsch parameterization is attainable locally in space 
[in the form (|60|)] [1C]. Additionally, an explicit construction of a, /?, and 9 can be given by 



the following pl| . 

Solve the equations 



which may also be presented as 



dx dy dz 

B X By B z 



£ ijk dr j B k = q ^ 62b ^ 

Solutions of these relations define two surfaces, called "magnetic surfaces", that are given by 
equations of the form 

S n (f) = constant (n = 1,2). (63) 

It follows from ( |62"| ) that these also satisfy 

B-V5 n = (n=l,2) (64) 

that is, the normals to S n are orthogonal to B, or B is parallel to the tangent of S n . The 
intersection of the two surfaces forms the so-called "magnetic lines", that is, loci that solve 
the dynamical system 

^l = B(r(r)) (65) 

where r is an evolution parameter. Finally, the Gaussian potentials a and (3 are constructed 
as functions of r only through a dependence on the magnetic surfaces, 

a{r) = a(S n (r)) 

(3{r) = P{S n (r)) (66) 



so that 

Vq x V/3 = (V5i x VS 2 ) £ ™^^. (67) 

Evidently as a consequence of (64), Va x V/? in (p^) is parallel to B, and because B is 
divergence-free a and f3 can be adjusted so that the norm of Va x V/3 coincides with \B\. 
Once a and (3 have been fixed in this way, 6 can easily be computed from A — a V/3. 

Neither the individual magnetic surfaces nor the Gauss potentials are unique. [By viewing 
A as a canonical 1-form, it is clear that the expression (|60[) retains its form after a canonical 
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transformation of a, (3.) One may therefore require that the Gaussian potentials a and (3 
simply coincide with the two magnetic surfaces: a = S\, (3 = S^- Nevertheless, for a given A 
and B it may not be possible to solve (]6^) explicitly. 

The Chern-Simons integrand A ■ B becomes in the Clebsch parameterization 

A - B = V6> • (Va x V/3) = V • {OB) = B-V9. (68) 

Thus having identified the Gauss potentials a and (3 with the two magnetic surfaces, we deduce 
from ( |64D and ( |68| ) three equations for the three functions (9, a, (3) that comprise the Clebsch 
parameterization. 

B ■ Va = B ■ V(3 = 

B ■ V# = Chern-Simons density A ■ B (69) 

Eq. (|68|) also shows that in the Clebsch parameterization the Chern-Simons density becomes 
a total derivative. 

A B = V (OB) (70) 

This does not mean that the Clebsch parameterization is unavailable when the Chern-Simons 
integral over all space is nonzero. Rather for a nonvanishing integral and well-behaved B 
field, one must conclude that the Clebsch function 9 is singular either in the finite volume of 
the integration region or on the surface at infinity bounding the integration domain. Then 
the Chern-Simons volume integral over (0) becomes a surface integral on the surfaces (d£l) 
bounding the singularities. 



drA-B= / dS-(0B) (71) 
n Jan 

Eq. (71) shows that the Chern-Simons integral measures the magnetic flux, modulated by 
and passing through the surfaces that surround the singularities of 9. 

The following explicit example illustrates the above points. 

Consider the vector potential whose spherical components are given by 

A r = (cos B)a'(r) 

= — (sin ©)— sin air) 
r 

A® = - (sine)- (l- cos a(r)). (72) 

(r, and 0, $ denote the conventional radial coordinate and the polar, azimuthal angles.) The 
function a{r) is taken to vanish at the origin, and to behave as 2ttv at infinity (y integer or 
half- integer). The corresponding magnetic field reads 

B r = -2(cos e)-j(l - cos o(r)) 

Bq = (sin 0)—a'(r) sin a(r) 
r 

Bq, = (sin e)-a'(r) (l - cos a(r)) (73) 
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and the Chern-Simons integral - also called the "magnetic helicity" in the electro dynamical 
context - is quantized (in multiples of 167r 2 ) by the behavior of a{r) at infinity 



f°° d 

drA-B = — 8-7T / dr — (air) — sin a(r)) 
Jo dr 

= -16ttV 



(74) 



In spite of the nonvanishing magnetic helicity, a Clebsch parameterization for ( |72| ) is readily 
constructed. In form notation, it reads 



(75) 



A = d(-2$)+2^1 - (sin 2 -) sin 2 0J d^ + tan" 1 (tan-) cos 9 
The magnetic surfaces can be taken from formula ( |75| ) to coincide with the Gauss potentials 

S\ = 2(1 — (sin 2 -) sin 2 0^ = constant 
S2 = & + tan" 1 (tan — ) cos 6 



constant 



(76) 



The magnetic lines are determined by the intersection of S\ and S^- 



cos — = ecos($ — ipo) 



sin 



1-e 2 



1 — e 2 cos 2 (<I> — ipo) 



(77) 



where e and tpo are constants. The potential 9 = — 2<3? is multivalued. Consequently the 
"surface" integral determining the Chern-Simons term reads 



dr A ■ B = / dr V • (-2&B) 



rdr / d6 5$ 



<E>=2tt 



(78) 



That is, the magnetic helicity is the flux of the toroidal magnetic field through the positive-x 
(x, z) half-plane. 



Problem 4 Consider a vector potential A, whose Clebsch parameterization reads Aj = di$ + 
cos 0<9j/i(r), where and $ are the azimuthal and polar angles of the vector r, and h is a 
nonsingular function of the magnitude of r. Show that the Chern-Simons density (magnetic 
helicity density) is given by AidjA^ = e l i k di&dj cos Qd^h^r). Consider the integral of the 
Chern-Simons density over all space. This integral may first be evaluated over a spherical ball 
0, and then the radius R of the ball is taken to infinity. When the integrand is a divergence 
of a vector, Gauss's theorem casts the volume integral onto a surface integral over the sphere 
dn bounding the ball: / = / n d 3 r V • V = J m dS -V = f^d<f> d© sin Gr 2 V r \ r=R , but 
singularities in V may modify the equality. The three derivatives in the above Chern-Simons 
density may be extracted in three different ways. Show that the result always is A7r[h(R) — 
h(0)] , but various singularities must be carefully handled. 
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Problem 5 Show that f d 3 r B ■ 5 A vanishes (apart from surface terms) where 5 A is a 
variation and A, B = V x A, as well as 5 A are presented in the Clebsch parameterization. 
When the variational principle is implemented by varying the components of A, one finds 
that I / d 3 r B 2 is stationary provided V x B = 0. Show that implementing the variation by 
varying the scalar functions in the Clebsch parameterization for A gives the weaker condition 
V x B = fj,B, where can depend on r. How is this r-dependence constrained? How is the 
constraint satisfied? 



There is another approach to the construction of (Abelian) vector potentials for which the 
(Abelian) Chern-Simons density is a total derivative, and as a consequence a Clebsch parame- 
terization for these potentials is readily found. The method relies on projecting an Abelian 
potential from a non- Abelian one, and it can be generalized to a construction of non- Abelian 
vectors for which the non-Abelian Chern-Simons density is again a total derivative. This will 
be useful for us when we come to discuss non-Abelian fluid mechanics. Therefore, I shall now 



explain this method - in its Abelian realization [12]. 

We consider an SU(2) group element g and a pure gauge SU(2) gauge field, whose matrix- 
valued 1-form is 



2i 

where a a are Pauli matrices. It is known that 



g -Ug = V a ?- (79) 



trGT 1 dg) 3 = -\ £ abc V a V b V c = -fyVV 3 (80) 
is a total derivative; indeed its spatial integral measures the winding number of the gauge 



function g [13]. Since V a is a pure gauge, we have 



dV a = -± £ abc V b V c (81) 

so that if we define an Abelian gauge potential A by projecting one SU(2) component of ([7S|) 
(say the third) A = V 3 , the Abelian Chern-Simons density for A is a total derivative, as is 
seen from the chain of equation that relies on (|80|) and (|8l|) . 



AdA = V 3 dV 3 = -V X V 2 V 3 = § ix(g~ x dg) 3 (82) 

Of course A = V 3 is not an Abelian pure gauge. 

Note that g depends on three arbitrary functions, the three SU(2) local gauge functions. 
Hence V 3 enjoys sufficient generality to represent the 3-dimensional vector A. Moreover, 
since A's Abelian Chern-Simons density is given by tr(g -1 dg) 3 , which is a total derivative, a 
Clebsch parameterization for A is easily constructed. We also observe that when the SU(2) 
group element g has nonvanishing winding number, the resultant Abelian vector possesses a 
nonvanishing Chern-Simons integral, that is, nonzero magnetic helicity. 
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Finally we remark that the example of a Clebsch-parameterized gauge potential A, pre- 
sented above in (|72[), is gotten by a projection onto the third isospin direction of a pure gauge 
SU(2) potential, constructed from group element g = exp((cr a /2z)f a a(r)) [12]. 

Further intricacies arise when the Clebsch parameterization is used in variational calcula- 
tions involving vector fields; see Ref. |ff3 . 
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3 Specific Models 

We now return to our irrotational models both relativistic and nonrelativistic, for which we 
shall specify an explicit force law and discuss further properties. 



3.1 Galileo-invariant nonrelativistic model 

Recall that the nonrelativistic Lagrangian for irrotational motion reads 

^Galileo = J • _ p (yfff _ y {p) ^ (g3) 

where V# = v. The Hamiltonian density TC is composed of the last two terms beyond the 
canonical 1-form J dr 9p, 

H = [ dr(p^P^ + V(p)) = [ drH. (84) 



2 

Various expressions for V appear in the literature. V(p) oc p n is a popular choice, appropriate 
for the adiabatic equation of state. We shall be specifically interested in the "Chaplygin gas" . 

V(p) = -, A > (85) 
P 

According to what we said before, the Chaplygin gas has enthalpy V = —X/p 2 , negative 
pressure P = —2X/p, and speed of sound s = V2X/p (hence A > 0). 

Chaplygin introduced his equation of state as a mathematical approximation to the physi- 
cally relevant adiabatic expressions with n > 0. (Constants are arranged so that the Chaplygin 
formula is tangent at one point to the adiabatic profile (llj.) Also it was realized that certain 
deformable solids can be described by the Chaplygin equation of state [16]. These days nega- 



tive pressure is recognized as a possible physical effect: exchange forces in atoms give rise to 
negative pressure; stripe states in the quantum Hall effect may be a consequence of negative 
pressure; the recently discovered cosmological constant may be exerting negative pressure on 
the cosmos, thereby accelerating expansion. 

For any form of V, the model possesses the Galileo symmetry appropriate to nonrelativistic 
dynamics. The Galileo transformations comprise the time and space translations, as well as 
space rotations. The corresponding constants of motion are the energy E. 

E = J drH = J dr{p^Y~ + V(p)j (time translation) (86) 

the momentum P (whose density "P equals the spatial current), 

P — j drV = j dr j = J dr pv (space translation) (87) 

and the angular momentum M !j . 

M ij = J dr (r^ - r 7 '^) (rotation) 
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The action of these transformations on the fields is straightforward: the time and space argu- 
ments are shifted or the space argument is rotated. Slightly less trivial is the action of Galileo 
boosts, which boost the spatial coordinate by a velocity u 

r -> R = r-tu. (89) 

The density field transforms trivially: its spatial argument is boosted, 

p(t,r)^p u (t,r)=p(t,R) (90) 

but the velocity potential acquires an inhomogeneous term. 

u 2 

0(t, r) -» 9 u (t, r) = 9(t, R) + u ■ r - —t (91) 

Those of you familiar with field theoretic realizations of the Galileo group will recognize the 
inhomogeneous term as the well-known Galileo 1-cocyle. It ensures that the velocity, the 
gradient of 9, transforms appropriately. 

v(t,r) ->v u (t, r) = v(t,r-tu) +u (92) 

The associated conserved quantity is the "boost generator". 

B = tP- j drpr (Galileo boost) (93) 

Finally, also conserved is the total number. 



N 



J dr p (particle number) (94) 



The corresponding transformation shifts 9 by a constant. 

These transformations and constants of motion fit into the general theory: the action is 
invariant against the transformations; Noether's theorem can be used to derive the constants 
of motion; their time independence is verified with the help of the equations of motion - 
indeed, the continuity equations (||), @, and (||) as well as the symmetry of T iJ guarantee 
this. Also, using the basic Poisson brackets (j25|) for the (p, 9) variables, one can check that each 
infinitesimal transformation is generated by Poisson bracketing with the appropriate constant; 
Poisson bracketing the constants with each other reproduces the Galileo Lie algebra with a 
central extension given by N, which corresponds to the familiar Galileo 2-cocycle. There are 
a total of \{d + + 2) Galileo generators in d space plus one time dimensions. Together 
with the central term, we have a total of \{d + + 2) + 1 generators. 

Another useful consequence of the symmetry transformations is that they map solutions of 
the equations of motion into new solutions. Of course, "new" solutions produced by Galileo 
transformations are trivially related to the old ones: they are simply shifted, boosted or rotated. 

[The free theory as well as the adiabatic theory with n = 1 + 2/d are also invariant 
against the SO(2,l) group of time translation, dilation, and conformal transformation Jl7|], 
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which together with the Galileo group form the "Schrodinger group" of nonrelativistic motion, 
whenever the energy-momentum "tensor" satisfies 2T°° = T l i [18].] 

But we shall now turn to the specific Chaplygin gas model, with V(p) = X/p, which 
possesses additional and unexpected symmetries. 

The Chaplygin gas action and consequent Bernoulli equation for the Chaplygin gas in (d, 1) 
space-time read 



I 



Chaplygin 



dt J dr(ep-\ P {ve) 



p 



X 



2 p 2 

This model possesses further space-time symmetries beyond those of the Galileo group Jl 
First of all, there is a one-parameter (u) time rescaling transformation 

t->T= e u t, 

under which the fields transform as 



(95) 
(96) 



0(t s r)-0 a ,(t,r) = e w 0(r j r) 
p(t,r)-> p u (t,r) = e- u p(T,r). 



(97) 

(98) 
(99) 



Second, in d spatial dimensions, there is a vectorial, (i-parameter space-time mixing trans- 
formation 



t -» T(t, r) = t + u-r + \oj 2 9{T, R 
r -» R{t,r) = r + u9(T,R) 



(100) 
(101) 



Note that the transformation law for the coordinates involves the 9 field itself. Under this 
transformation, the fields transform according to 



6(t,r) -> u (t,r) = 6(T,R) 
p(t,r)^p„(t,r) = p(T,R) 



\J\' 



(102) 
(103) 



J = det 



l-u-Ve(T,R)-—9(T,R) 



(104) 



with J the Jacobian of the transformation linking (T, R) — > (t, r) 

( dT dT\ 

dt dri 
OR 1 dRj- 
\~dt ~dri ) 

(The time and space derivatives in the last element are with respect to t and r.) One can again 
tell the complete story for these transformations: The action is invariant; Noether's theorem 
gives the conserved quantities, which for the time rescaling is 



D = tH 



dr p9 (time rescaling) 



(105) 
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while for the space-time mixing one finds 

G = j dr (rH — 9V) (space-time mixing). (106) 

The time independence of D and G can be verified with the help of the equations of motion 
(continuity and Bernoulli). Poisson bracketing the fields 9 and p with D and G generates the 
appropriate infinitesimal transformation on the fields. 

So now the total number of generators is the sum of the previous (d + 1) (d + 2) + 1 with 
1 + d additional ones 

\{d + l)(d + 2) + 1 + 1 + d = + 2)(d + 3). (107) 

When one computes the Poisson brackets of all these with each other one finds the Poincare 
Lie algebra in one higher spatial dimension, that is, in 1) -dimensional space-time, where 

the Poincare group possesses ^(d+2)(d+3) generators |2C|| . Moreover, one verifies that (t, 9, r) 
transform linearly as a (d + 2) Lorentz vector in light-cone components, with t being the "+" 
component and 6 the "— " component [21|. 

Presently, we shall use these additional symmetries to generate new solutions from old ones, 
but, in contrast with what we saw earlier, the new solutions will be nontrivially linked to the 
former ones. Note that the additional symmetry holds even in the free theory. 

Before proceeding, let us observe that p may be eliminated by using the Bernoulli equation 
to express it in terms of 9. In this way, one is led to the following p- independent action for 9 
in the Chaplygin gas problem: 



jChaplygin = _ 2 ^ J dt J dr ^9 + ^31, (1 08 ) 

Although this operation is possible only in the interacting case, the interaction strength is seen 
to disappear from the equations of motion. 

8 1 h V , V " (109) 



A merely serves as an overall factor in the action. 

The action ( JlOSj ) looks unfamiliar; yet it is Galileo invariant. [The combination 9 + ^(V#) 2 
responds to Galileo transformations without a 1-cocycle; see (|9l|).1 Also /^ haplygm possesses 
the additional symmetries described above, with 9 transforming according to the previously 
recorded equations. 

Let us discuss some solutions. For example, the free theory is solved by 

2 

9 M = Y t ( n °) 



which corresponds to the velocity 



v(t,r) = ~. (Ill) 



28 
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Galileo transforms generalize this in an obvious manner into a set of solutions. (The charge 
density p is determined by its initial condition. In the free theory, p is an independent quantity, 
and I shall not discuss it here.) Performing on the above formula for 9 the new transformations 
of time-rescaling and space-time mixing, we find that the solution is invariant. 

We can find a solution similar to ( |110| ) in the interacting case, for d > 1, which we henceforth 
assume (the d = 1 case will be separately discussed later). One verifies that a solution is 

s(4 '-» = -wrva *-> = \/f( d - 1 )7 < 112 > 

v(t,r) =- —— j(t,r) = -s(t)^—f. (113) 
Note that the speed of sound 

V2X ^fdr 



p (a — l)t 

exceeds v. Again this solution can be translated, rotated, and boosted. Moreover, the solution 
is time-rescaling-invariant. However, the space-time mixing transformation produces a wholly 
different kind of solution. This is best shown graphically, where the d = 2 case is exhibited 



(see figure) [22] 



Another interesting solution, which is essentially one-dimensional (lineal), even though it 
exists in arbitrary spatial dimension, is given by 

0(t,r) = e(h-r)+u-r-±t(u 2 - {h-uf). (115) 

Here h is a spatial unit vector, u is an arbitrary vector with dimension of velocity, while G is an 
arbitrary function with static argument, which can be boosted by the Galileo transform (|9l|). 
The corresponding charge density is time-independent. 



p(t, r) = (116) 

rv ' n-u + Q'(n • r) v ' 

The current is static and divergenceless. 

i— ( u — h(h ■ u) \ 

j(t,r) = V2\(h+- K > 117 

\ n ■ r + G'(n • r)J 

The sound speed s = V2X/p = Q'(n ■ r) + n ■ u is just the h component of the velocity 
v = WO = nG'(n • r) + u. 

Finally, we record a planar static solution to ( |109| ) 

0(r) = 8(ni -r/n 2 • r) (118) 

where h\ and fi2 are two orthogonal unit vectors |23||. 



Problem 6 Show that the solution for 6 given in ( |110[) is invariant under the time rescaling 
transformation (|97|), (p8|), and under the space-time mixing transformation ( |10C| )-( |lO^ ). 
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10 10 

The transformed density p(t,r). 
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3.2 Lorentz-invariant relativistic model 

We now turn to a Lorentz-invariant generalization of our Galileo-invariant Chaplygin model 
in (d, 1) -dimensional space-time. We already know from (|34|)-(|37|) how to construct the free 
Lagrangian by using Eckart's method with a relativistic kinetic energy. 



T (v) = -cVl - v 2 /c 2 (119) 

Recall that mass has been scaled to unity, and that we retain the velocity of light c to keep 
track of the nonrelativistic c —* oo limit. Evidently, the momentum is 

dT(v) v , . 

P=^r —= - 120) 

Thus the free relativistic Lagrangian, with current conservation enforced by the Lagrange 
multiplier 6, reads [compare 



Z Lorentz = J dr (- C 2 p^ 1 - V 2 / C 2 + 6 (p + V ■ (pv))^. (121) 

This may be presented in a Lorentz-covariant form in terms of a current four-vector = 
(cp,pv). LQ° rentz of equation ( |121| ) is thus equivalent to [compare (]56|), (|57|)] 



^Lorcntz = J dr (_^ dti _ cy^). (122) 



Eliminating v in ( |l2l| ), we find, as before, that 



P = d ^= - V =V9, v= - V ° , (123) 
dv Vl^/c 2 y/l + (V9) 2 /c 2 K ' 



and the free Lorentz-invariant Lagrangian reads [compare (|3q), (p7|)1 

Ll orcntz = j dr(6p- pc 2 y/TTWWJ^)- (124) 

To find L[f alileo of @) as the nonrelativistic limit of £L° rcntz in (jl2|), a nonrelativistic 6 
variable must be extracted from its relativistic counterpart. Calling the former #nr and the 
latter, which occurs in (124), we define 

R = -c 2 t + # N R. (125) 

It then follows that apart from a total time derivative 

L Lorentz > L? alUe °. (126) 

c^oo 

Next, one wants to include interactions. While there are many ways to allow for Lorentz- 
invariant interactions, we seek an expression that reduces to the Chaplygin gas in the nonrel- 
ativistic limit. Thus, we choose [|4j] 

L^^j^-^^T^^TWW), (127) 
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where a is the interaction strength. [The reason for the nomenclature will emerge presently.] 
We see that, as c — » oo, 

^Born-Infeld > ^Chaplygin (128) 

a c^oo A=a 2 /2 ' ^ > 

[Again is extracted from #r as in (|125| ) and a total time derivative is ignored.] Although 



it perhaps is not obvious, (127) defines a Poincare-invariant theory, and this will be explicitly 
demonstrated below. Therefore, ^Born-infeid p OSSesses Lorentz and Poincare symmetries in 
(d, 1) space-time, with a total of | (d + 1) (d + 2) + 1 generators, where the last " + 1" refers to 
the total number N = J dr p. 

When a = 0, the model is free and elementary. It was demonstrated previously [eqs. (ET~ 



(42)] that the free equations of motion are precisely the same as in the nonrelativistic free 
model, so the complete solution ([43|)-(|4"7"|) works here as well. For a 7^ 0, in the presence of 
interactions, one can eliminate p as before, and one is left with a Lagrangian just for the 6 
field. It reads 

L Born-Infcld = _ a J dr ^2 _ ^)2. ( 129 ) 

This is a Born-Infeld-type theory for a scalar field 6; its Poincare invariance is manifest, and 
again, the elimination of p is only possible with nonvanishing a, which however disappears 
from the dynamics, serving merely to normalize the Lagrangian. 



The equations of motion that follow from ( 127 ) read 



' +v ' W^fW =° (130> 



s + ,Af±(Zf = „. (131) 
V p c +a 



The density p can be evaluated in terms of 9 from (|131|); then Q13C|) reads 



d a [ - ^d a 6)=0 (132) 

V c2 -(W 



which also follows from (129). After is extracted from #r as in ( |125|) we see that in the 



nonrelativistic limit ^Bom-Meld Q Q r O becomes L^hapiygin q{ q Qr 



rBorn-Infeld t Chaplygin /ino\ 
L \=cfll2 \ 166 

c — >oo /* 



and the equations of motion ( |130| )-( ^32] ) reduce to ([T|), (|96[), and ( 109 ). 



In view of all the similarities to the nonrelativistic Chaplygin gas, it comes as no surprise 
that the relativistic Born-Infeld theory possesses additional symmetries. These additional 
symmetry transformations, which leave (|127| ) or (|l29j ) invariant, involve a one-parameter (uj) 
reparameterization of time, and a (i-parameter (u>) vectorial reparameterization of space. Both 
transformations are field dependent. 
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The time transformation is given by an implicit formula involving also the field 9 [25] 

t 6(T,r) 



t->T(*,r) = 
while the field transforms according to 

0(t,r)-0 w (t,r) 



cosh c 2 uj 



+ 



9(T,r) 



tanh c 2 u> 



c 2 itanhc 2 u; . 



(134) 



(135) 



cosh C 2 LO 

[We record here only the transformation on 9; how p transforms can be determined from the 
(relativistic) Bernoulli equation, obtained by varying p in ( |127| ), which expresses p in terms 
of 9. Moreover, (|135| ) is sufficient for discussing the invariance of ( |129|) ,1 The infinitesimal 



generator, which is time independent by virtue of the equations of motion, is [26] 
D = f dr (cHp + 9^fp 2 c 2 + a 2 ^c 2 + {V9) 2 ) 

A second class of transformations involving a reparameterization of the spatial variables is 



implicitly defined by [25|. 



„ / x „ „ / „n tan cm . , A 
r -> R(t, r) = r - u9(t, R) h w(c2) • r) 



1 — COS CUJ 



COS CCJ 



(137) 



„ , . . . . 9{t, R) — c(uj ■ r) sinew „ nnl 
9(t, r) -» w (i, r) = v y i ^ (138) 

COS CUJ 

Cj is the unit vector u/uj and u = ■ u. The time-independent generator of the infinitesimal 
transformation reads [26] 

G = j dr(c 2 rp + 9 P V9) 

= / ir(c'r P + m (space re P ara m e,en Z a t i„„). (139) 



Of course the Born-Infeld action ( |127| ) or ( |12S| ) is invariant against these transformations, 
whose infinitesimal form is generated by the constants. 

With the addition of D and G to the previous generators, the Poincare algebra in (d+ 1, 1) 
dimension is reconstructed, and (t,r,9) transforms linearly as a (d + 2)-dimensional Lorentz 
vector (in Cartesian components) [ffif| . Note that this symmetry also holds in the free, a = 0, 
theory. 

It is easy to exhibit solutions of the relativistic equation (132), which reduce to solutions of 
the nonrelativistic, Chaplygin gas equation ( 109[ ) [after — c 2 t has been removed, as in ( |125| )]. 
For example 



9(t,r) = -c\/cH 2 + — - (140) 
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solves (132) and reduces to (112). The relativistic analog of the lineal solution (115) is 



9(t, r) = 0(n • r) + n • r — ct\J c 2 + u 2 — (n ■ u) 2 . 



(141) 



[Note that the above profiles continue to solve ( |132j ), even when the sign of the square root is 
reversed, but then they no longer possess a nonrelativistic limit.] 

Additionally there exists an essentially relativistic solution, describing massless propaga- 
tion in one direction: according to (132), 9 can satisfy the wave equation D9 = 0, provided 
(d^O) 2 = constant, as for example with plane waves 



0(t,r) = f(n-r±ct) 



(142) 



where (8^0) 2 vanishes. Then p reads, from (131), 



P = T-nf 



(143) 



Other solutions are given in ref. [27]. 



3.3 Some remarks on relativistic fluid mechanics 

The Born-Infeld model reduces in the nonrelativistic limit to the Chaplygin gas. Equations 
governing the latter belong to fluid mechanics, but the Born-Infeld equations do not readily 
expose their fluid mechanical structure. Nevertheless they do in fact describe a relativistic 
fluid. In order to demonstrate this, we give a precis of relativistic fluid mechanics. 

Usually the dynamics of a relativistic fluid is presented in terms of the energy-momentum 
tensor, 9^ u , and the equations of motion are just the conservation equations d^9^ v = p§|]. 
[We denote the relativistic energy- momentum tensor by 9^ u , to distinguish it from the nonrel- 
ativistic T^ v introduced in (Q) and (||). The limiting relation between the two is given below.] 
But we prefer to begin with a Lagrange density 



£ = -i\-/(VHi)- ( 144 ) 

Here is the current Lorentz vector = (cp,j). The comprise a set of auxilliary variables; 
in the relativistic analog of irrotational fluids we take = 8^9, more generally 

a M = 8^9 + ad„/3 (145) 



so that the Chern-Simons density of ctj is a total derivative [compare (|56[), (57)]. The function / 
depends on the Lorentz invariant j'^j^ = (?p 2 —j 2 and encodes the specific dyanamics (equation 
of state). 

The energy momentum tensor for C is 

V = ~9^C + Jf^f\^pr a ). (146) 



34 
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[One way to derive (146) from ( |144| ) is to embed that expression in an external metric tensor 
gaui which is then varied; in the variation and are taken to be metric-independent and 
jfj, = 9iivj v ■] Furthermore, varying in ( |144| ) shows that 

a^-^^nVFL) (147) 

so that (146) becomes 

V = -9^[nf[n) - f{n)} + u^u u nf\n) (148) 

where we have introduced the proper velocity Un by 

= nu^ u% = 1 (149) 

so that n is proportional to the proper density and 1/n is proportional to the specific volume. 
Eq. (148) identifies the proper energy density e and the pressure P (which coincides with £) 
through the conventional formula |2§|]. 



O llv = -g llv P + u ll u v (P + e) (150) 

Therefore, in our case 

e = f(n) (151) 

P = nf{n) - f(n) . (152) 

The thermodynamic relation involving entropy S reads 

pd Q +d G) Kds < i53 > 

where the proportionality constant is determined by the temperature. With ( |151|) and (|152j ) 
the left side of (|153| ) vanishes and we verify that entropy is constant, that is, we are dealing 
with an isentropic system, as has been stated in the very beginning. 
For the free system, the pressure vanishes, so we choose f(n) = cn 

Co = -fa^ - Cy qj%. (154) 

In the "irrotational" case, a M = 8^8, and with j = pv this Lagrange density produces the free 
^Lorentz of ( gg^ Q ( apart f rom a total time derivative). 

For the Born-Infeld theory, we present the pressure P in ( |152|) by choosing f{n) = cV a 2 + n 2 , 
which corresponds to the pressure P = —a 2 c/\/a 2 + n 2 . When the current is written as 

lu = — = 155a) 

^ ^c 2 - (pp 



so that 



n = a \l c 2- W ) 2 (155b) 
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the Lagrange density 



a?c 



coincides with that for the Born-Infeld model in eq. ( |129| ). 

Other forms for / give rise to relativistic fluid mechanics with other equations of state. 

It is interesting to see how the nonrelativistic limit of 9^ v in ( |148[) produces T^ u of (f|)-(||). 
It is especially intriguing to notice that 9^ v is symmetric but T^ v is not. To make the connec- 
tion we recall that = 1/^/1 — v 2 /c 2 (l, v/c), we observe that n = \J p 2 c 2 — j 2 , set j = pv and 
conclude that n = 



pc^J\ — v 2 /c 2 ~ pc — {pv 2 /2c). Also f(n) is chosen to be cn + V(n/c), and thus P = 
nf{n) - f{n) = (n/c)V'(n/c) - V(n/c). It follows that 

nc - (v 2 n/c 3 )V pc 2 - pv 2 /2 

V = i 27~2 r V "1 2T^~ + v W 

1 — v z /c z 1 — v z /c z 

2 

~PC 2 + ^~ + V(p) = pc 2 + T°° . (156) 

Thus, apart from the relativistic "rest energy" pc 2 , 9°° passes to T°°. The relativistic energy 
flux is dp (because ^9^° = \9°° + d^°) 

n T //\ 7 P° 2 ~ P v2 / 2 + PV'(P) 



c9 3 ° = - „[nc+-V 



1 — v 2 /c 2 v c / 1 — v 2 /c 2 

j J c 2 + /V (u 2 /2 + V'(p)) = j j c + r^'° . (157) 

Again, apart from the 0(c 2 ) current, associated with the 0(c 2 ) rest energy in 9°°, is 
obtained in the limit. The momentum density is 9 m /c (because 9^ u has dimension of energy 
density). Thus 

9 oi /c = - Vt /f. _ (nc + -V') sa /w* = V i . (158) 



1 — v 2 /c 2 V c 
Finally, the momentum flux is obtained directly from 9 lJ . 



0* = S nly'_y\ + ( nc+ _y> 

\ c / c z — 7J Z V c / 

» ^' (pV'(p) - V(/j)) + uV> = r ij (159) 



From the limiting formula n ~ /;c we also see that the pressure in ( 1553 ) tends to the 
Chaplygin expression —a/ p. 
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4 Common Ancestry: The Nambu-Goto Action 



The "hidden" symmetries and the associated transformation laws for the Chaplygin and Born- 
Infeld models may be given a coherent setting by considering the Nambu-Goto action for a 
d-brane in (d + 1) spatial dimensions, moving on (d + 1, l)-dimensional space-time. In our 
context, a d-brane is simply a ci-dimensional extended object: a 1-brane is a string, a 2-brane 
is a membrane and so on. A d-brane in (d + 1) space divides that space in two. 
The Nambu-Goto action reads 



Here X^ is a (d + 1, 1) target space-time (d-brane) variable, with fi extending over the range 
/x = 0, 1, . . . , d, d + 1. The ip a are "world-volume" variables describing the extended object 
with a ranging a = 0, 1, . . . , d; ip a , a = 1, . . . , d, parameterizes the (i-dimensional d-brane that 
evolves in ip°. 

The Nambu-Goto action is parameterization invariant, and we shall show that two different 
choices of parameterization ( "light-cone" and "Cartesian" ) lead to the Chaplygin gas and Born- 
Infeld actions, respectively. For both parameterizations we choose (X 1 , . . . , X d ) to coincide 
with (if 1 , . . . , tp d ), renaming them as r (a ci-dimensional vector). This is usually called the 
"static parameterization" . (The ability to carry out this parameterization globally presupposes 
that the extended object is topologically trivial; in the contrary situation, singularities will 
appear, which are spurious in the sense that they disappear in different parameterizations, and 
parameterization-invariant quantities are singularity- free.) 

4.1 Light-cone parameterization 

For the light-cone parameterization we define X ± as ^(X° ± X d+1 ). X + is renamed t and 
identified with y/2\ ip° . This completes the fixing of the parameterization and the remaining 
variable is X~ , which is a function of y? and ip, or after redefinitions, of t and r. X~ is 
renamed as 0(t,r) and then the Nambu-Goto action in this parameterization coincides with 
the Chaplygin gas action I^p^ 11 [ n flfij|) 

4.2 Cartesian parameterization 

For the second, Cartesian parameterization X° is renamed ct and identified with ap . The re- 
maining target space variable X d+1 , a function of </?° and ip, equivalently of t and r, is renamed 
9(t,r)/c. Then the Nambu-Goto action reduces to the Born-Infeld action f dt L^ orn " Infeld , 




(160) 



G = (-l) d det 



dXt 1 dXjt 
dip a dpi 3 



(161) 
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4.3 Hodographic transformation 

There is another derivation of the Chaplygin gas from the Nambu-Goto action that makes use of 
a hodographic transformation, in which independent and dependent variables are interchanged. 
Although the derivation is more involved than the light-cone/static parameterization used in 
Section |4.1| above, the hodographic approach is instructive in that it gives a natural definition 
for the density p, which in the above static parameterization approach is determined from 9 
by the Bernoulli equation (p6|). 

We again use light-cone combinations: -^(X° + X d+1 ) is called r and is identified with 
ip°, while -^(X° — X d+1 ) is renamed 6. At this stage the dependent, target-space variables 
are 6 and the transverse coordinates X : X 1 , (i = 1, . . . , d), and all are functions of the world- 
volume parameters ip° = T and if: ip r , (r = 1, . . . ,d);d T indicates differentiation with respect 
to r — (p® , while d T denotes derivatives with respect to (p r . The induced metric G a p — Sy" dipP 
takes the form 

G = (g 00 g os \ = ( 2d T e - {d T xf d s e - d T x ■ d s x\ 

aP \G ro -g TS ) \d T e-d r X-d T X -d r X ■ d s X J '~ ! 

The Nambu-Goto Lagrangian now leads to the canonical momenta 



p (163a) 
n (163b) 



dd T X 

OCng 
dd T 9 

and can be presented in first-order form as 

£ng = P ■ d T X + Ud T + ^(p 2 + g)+ u r (p ■ d r X + Ild r 6) (164) 
where g = det g rs and 

u r = d T X ■ d r X - d r 6 (165) 
acts as a Lagrange multiplier. Evidently the equations of motion are 

d T X = -ip - u r d r X (166a) 

d T 9 = ^(p 2 + g)-u r d r 8 (166b) 

d T p = ~d r ^gg rs d s X^j - d r (u r p) (166c) 
d T U = -d r (u r U) (166d) 

Also there is the constraint 



p ■ d r x + m r e = o 



(167) 
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[That u r is still given by (165) is a consequence of ( |166a| ) and (167).] Here g rs is inverse 



to g rs , and the two metrics are used to move the (r, s) indices. The theory still possesses 
an invariance against redefining the spatial parameters with a r-dependent function of the 
parameters; infinitesimally: 5ip r = — f r (r, cp), 59 = f r d r 9, 5X l = f r d r X l . This freedom may 
be used to set u r to zero and II to —1. 

Next the hodographic transformation is performed: Rather than viewing the dependent 
variables p, 9, and X as functions of r and <p, X(r,(p) is inverted so that (p becomes a 
function of t and X (renamed t and r, respectively), and p and 9 also become functions of t 
and r. It then follows from the chain rule that the constraint (|167|) (at II = —1) becomes 



= S^"^) (168) 



dip r \ dX 
and is solved by 

p = V9 . (169) 

Moreover, according to the chain rule and the implicit function theorem, the partial derivative 
with respect to r at fixed ip [this derivative is present in ( |164| )1 is related to the partial derivative 
with respect to r at fixed X = r by 

d T = ^ + V9 • V (170) 

where we have used the new name "i" on the right. Thus the Nambu-Goto Lagrangian - the 
(p integral of the Lagrange density (|164| ) (at u r = 0, II = — 1) - reads 

L NG = J dip{p -V9-9-V9-V9- \{p 2 + g)}. (171a) 



But use of ( |169[) and of the Jacobian relation dip = dr det = shows that 



L NG = J H-^Q ~ ^ W ) 2 " (171b) 

With the definition 

^g-=V2\/p (171c) 
Lnc7 becomes, apart from a total time derivative 

L * G = 7Txf dr { P ~ ^ (V0)2 ~ (171d) 
Up to an overall factor, this is just the Chaplygin gas Lagrangian in (|95|). 

The present derivation has the advantage of relating the density p to the Jacobian of 
the X 1 — > ip transformation: p = V2Xdet^i. (This in turn shows that the hodographic 
transformation is just exactly the passage from Lagrangian to Eulerian fluid variables - a 
remark aimed at those who are acquainted with the Lagrange formulation of fluid motion ||.) 

Finally, let me observe that the expansion of the Galileo symmetry in (d, 1) space-time to 
a Poincare symmetry in (d + 1, 1) space-time can be understood from a Kaluza-Klein-type 
framework [30]. 
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4.4 Interrelations 

The relation to the Nambu-Goto action explains the origin of the hidden (d + 1, 1) Poincare 
group in our two nonlinear models on (d, 1) space-time: Poincare invariance is what remains 
of the reparameterization invariance of the Nambu-Goto action after choosing either the light- 
cone or Cartesian parameterizations. Also the nonlinear, field dependent form of the transfor- 
mation laws leading to the additional symmetries is understood: it arises from the identification 
of some of the dependent variables (X^) with the independent variables {(p a ). 

The complete integrability of the d = 1 Chaplygin gas and Born-Infeld model is a con- 
sequence of the fact that both descend from a string in 2-space; the associated Nambu-Goto 
theory being completely integrable. We shall discuss this in Section ||. 

We observe that in addition to the nonrelativistic descent from the Born-Infeld theory to 
the Chaplygin gas, there exists a mapping of one system on another, and between solutions 
of one system and the other, because both have the same d-brane ancestor. The mapping 
is achieved by passing from the light-cone parameterization to the Cartesian, or vice-versa. 
Specifically this is accomplished as follows: 

Chaplygin gas — > Born-Infeld: Given djs? R (t,r), a nonrelativistic solution, determine 
T(t, r) from the equation 



T+\e NR (T,r) = V2t (172) 

Then the relativistic solution is 



6 R (t, r) = -±=c 2 T - -j=0 N r(T, r) = c 2 (V2T - t) (173) 

Born-Infeld — > Chaplygin gas: Given 6 R (t,r), a relativistic solution, find T(t,r) from 

T + \e R (T,r) = V2t (174) 



Then the nonrelativistic solution is 



9 NR (t, r) = -Lc 2 T - ^=0r(T, r) = c 2 (V2T - t) (175) 

The relation between the different models is depicted in the figure below. 

As a final comment, I recall that the elimination of p, both in the nonrelativistic (Chap- 
lygin) and relativistic (Born-Infeld) models is possible only in the presence of interactions. 
Nevertheless, the ^-dependent (p-independent) resultant Lagrangians contain the interaction 
strengths only as overall factors; see (108) and ( |129| ). It is these ^-dependent Lagrangians that 



correspond to the Nambu-Goto action in various parameterizations. Let us further recall the 
the Nambu-Goto action also carries an overall multiplicative factor: the d-brane "tension", 
which has been suppressed in ( |160| ). Correspondingly, for a "tensionless" d-brane, the Nambu- 
Goto expression vanishes, and cannot generate dynamics. This suggests that an action for 
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Dualities and other relations between nonlinear equations. 



"tensionless" d-branes could be the noninteracting fluid mechanical expressions fl95[) , (|127| ), 
with vanishing coupling strengths A and a, respectively. Furthermore, we recall that the non- 
interacting models retain the higher, dynamical symmetries, appropriate to a d-brane in one 
higher dimension. 
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5 Supersymmetric Generalization 

Once proven the fact that a bosonic Nambu-Goto theory gives rise to and links together the 
Chaplygin gas and Born-Infeld models, which are irrotational in that the velocity of the former 
and the momentum of the latter are given by a gradient of a potential, one can ask whether 
there is a d-brane that produces a fluid model with nonvanishing vorticity. 

We shall show that this indeed can be achieved if one starts with a super-d-brane, and 
moreover the resulting fluid model possesses super symmetry. However, since theories of ex- 
tended "super" objects cannot be formulated in arbitrary dimensions, we shall consider the 
fluid in two spatial dimensions, namely, on the plane plj . 



5.1 Chaplygin gas with Grassmann variables 



We begin by positing the fluid model. The Chaplygin gas Lagrangian in (95) is supplemented 
by Grassmann variables ip a that are Majorana spinors [real, two-component: i/j* = ip a , a = 1,2, 

(VlV>2)*=^]- 

The associated Lagrange density reads 

\ fey \ 

C = -p(9 - - 3P(V0 - ^V>VV0 2 — ipa ■ VV> . (176) 

Here a 1 are two (i = 1,2), 2 x 2, real symmetric Dirac "alpha" matrices; in terms of Pauli 
matrices we can take a 1 = a 1 , c? = a 3 . Note that the matrices satisfy the following relations, 
which are needed to verify subsequent formulas 

a ab a bc = SijS ac ~ £ ij £ac 

a ab a ld = $ac5bd - &ab&cd + Sad^bc, (177) 

e a b is the 2x2 antisymmetric matrix e = ia 2 . In equation ( |176[) A is a coupling strength which 
is assumed to be positive. The Grassmann term enters with coupling V2X, which is correlated 
with the strength of the Chaplygin potential V(p) = X/p in order to ensure super symmetry, 
as we shall show below. It is evident that the velocity should be defined as 

v = ve- ±v>vv> . (178) 

The Grassmann variables directly give rise to a Clebsch formula for v, and provide the 
Gauss potentials. The two-dimensional vorticity reads u = e l ^diiP = —he^di^djifj = — \ VV> X 



V^. The variables {9, p} remain a canonical pair, while the canonical 1-form in ( 176 ) indicates 
that the canonically independent Grassmann variables are y/pip so that the antibracket of the 

ip's is 

{Mr),Mr')} = -^f L J(r-r') . (179) 
p(r) 
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One verifies that the algebra (|4|) or ( |i~6| ) is satisfied, and further, one has 

{9{r)^(r)} = -^—^ r )6(r - r') (180) 

{v(r), V>(r')} = -^P^5(r - r') (181) 
p(r) 

{V(r),ip(r')} = -Vip(r)5(r - r') . (182) 

The momentum density "P is given by the bosonic formula "P = pv, but the Grassmann 
variables are hidden in v, by virtue of (|178j ). 
The equations of motion read 



p + V-(pv) = (183) 
7 + ^P 



+ v ■ V6 = \v 2 + — + ipa ■ Vt/> (184) 



v2A 

i> + v ■ W = - — « ■ (185) 
P 

and together with ( |1T8|) they imply 

A v^2~A 

i, + u . = + - — (Vip)a-Vip. (186) 
All these equations may be obtained by bracketing with the Hamiltonian 

H = J d 2 r^pv 2 + ^ + ^^^a- w) = Jd 2 rft (187) 

when (|D, (ID as well as (|l79|)-(|l8ll) are used. 

We record the components of the energy-momentum "tensor" , and the continuity equations 
they satisfy. The energy density £ = T°°, given by 

£ = \pv 2 + - + ^ ^ a • VV = T°° (188) 
p A 

satisfies a continuity equation with the energy flux T 3 °. 



T jo = pv j ( \v 2 - 4) + — ■ - -ipdjib + -e jk ^ed k ^ (189) 
vz p 1 ' 2 p J p 

f°° + djT jo = (190) 

This ensures that the total energy, that is, the Hamiltonian, is time-independent. Conservation 
of the total momentum 

d 2 rV (191) 
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follows from the continuity equation satisfied by the momentum density V 1 
momentum flux, that is, the stress tensor T lJ . 



pv l v 3 



v p 



, _ A v2A , • 
- — tpa ■ vipj H — r ipa J Oi'ip 



f oi + djT ji = 



T oi and the 

(192) 
(193) 



But T %3 is not symmetric in its spatial indices, owing to the presence of spin in the problem. 
However, rotational symmetry makes it possible to effect an "improvement", which modifies 
the momentum density by a total derivative term, leaving the integrated total momentum 
unchanged (provided surface terms can be ignored) and rendering the stress tensor symmetric. 
The improved quantities are 



V\ = Tf = pv i + \e i id j {pipei>) 



T i3 



pv l v 3 



2A 



ipot ■ Vf/j ) + 



v p 

- \d k \(e ki v j + e^v^pnpeip 



(194) 

(195) 
(196) 



Tf + djT/ = . 

It immediately follows from the symmetry of Tj J that the angular momentum 

M = j d 2 re lj r i Vj = J d 2 r pe ij r i v j + \ J d 2 r pipeip 

is conserved. The first term is clearly the orbital part (which still receives a Grassmann 
contribution through v), whereas the second, coming from the improvement, is the spin part. 



(197) 



Indeed, since |e 



E, we recognize this as the spin matrix in (2+1) dimensions. 
The extra term in the improved momentum density ( |194| ), ^e i3 dj(ptp£tp), can then be readily 
interpreted as an additional localized momentum density, generated by the nonhomogeneity of 
the spin density. This is analogous to the magnetostatics formula giving the localized current 
density j m in a magnet in terms of its magnetization m: j m = V x m. All in all, we are 
describing a fluid with spin. 
Also the total number 



N 



d 2 r p 



(198) 



is conserved by virtue of the continuity equation (|183j ) satisfied by p. Finally, the theory is 
Galileo invariant, as is seen from the conservation of the Galileo boost, 



B = tP 



d 2 r rp 



(199) 



which follows from (|183| ) and ( 191 ). The generators H, P, M, B and close on the (extended) 
Galileo group. [The theory is not Lorentz invariant in (2 + l)-dimensional space-time, hence 
the energy flux T JO does not coincide with the momentum density, improved or not.] 
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We observe that p can be eliminated from ( |176| ) so that C involves only 6 and ip. From 
( gH) and (pi) it follows that 



p = yf\(e - + \v 2 ) 1/2 . (200) 

Substituting into ( |176| ) produces the supersymmetric generalization of the Chaplygin gas La- 
grange density in (|108|). 



c = -2^^ 26 -ipip + (ye - i^V^) 2 + lipa ■ VV>} (201) 

Note that the coupling strength has disappeared from the dynamical equations, remaining 
only as a normalization factor for the Lagrangian. Consequently the above elimination of p 
cannot be carried out in the free case, A = 0. 

5.2 Supersymmetry 

As we said earlier, this theory possesses supersymmetry. This can be established, first of all, 
by verifying that the following two-component supercharges are time-independent Grassmann 
quantities. 



Qa = j d 2 r [pw • (a ab i{j b ) + V2X^a] ■ (202) 



Taking a time derivative and using the evolution equations ( |183| )-( |i"86| ) establishes that Q a = 0. 

Next, the supersymmetric transformation rule for the dynamical variables is found by 
constructing a bosonic symmetry generator Q, obtained by contracting the Grassmann charge 
with a constant Grassmann parameter rj a , Q = ij a Q a , and commuting with the dynamical 
variables. Using the canonical brackets one verifies the following field transformation rules: 

5p = {Q,p} = -V-p(r ] arP) (203) 



56 = {Q, 6} = -foaif)) ■ V6 - j(ryaV) • ipVip + rtf (204) 



2A 
2j 

a/2A 

8tJ) = {Q, ip} = —(ijatf;) ■ Vt/> — v ■ our] r\ (205) 



v 2A 

5v = {Q,v} = -(rjaip) ■ Vd + — —r/Vip . (206) 

Supersymmetry is reestablished by determining the variation of the action J dtd 2 rC conse- 
quent to the above field variations: the action is invariant. One then reconstructs the super- 



charges ( 202 ) by Noether's theorem. Finally, upon computing the bracket of two supercharges, 
one finds 



MQa,v b 2Qb} = 2(mm)H 



(207) 
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which again confirms that the charges are time-independent: 

{H,Q a } = 0. (208) 

Additionally a further, kinematical, supersymmetry can be identified. According to the 
equations of motion the following two supercharges are also time-independent: 

Qa = J d 2 rp^ a . (209) 

Q = fj a Q a effects a shift of the Grassmann field: 

Sp = {Q,p} = (210) 

60 = {Q,9} = -§(#) (211) 

51> = {Q,M = -fj (212) 

Sv = {Q,v} = . (213) 



This transformation leaves the Lagrangian invariant, and Noether's theorem reproduces (|209| ) . 
The algebra of these charges closes on the total number N. 

{fftQa,v b 2 Q b } = (mm)N (214) 

while the algebra with the generators ( |202| ) , closes on the total momentum, together with a 
central extension, proportional to volume of space $7 = J d 2 r 

{tQa, V b Qb} = (TO) -P + V2X {f}er])n . (215) 

The supercharges Q a ,Qa, together with the Galileo generators (H, P, M, and B), with N 
form a superextended Galileo algebra. The additional, nonvanishing brackets are 

{M,Q a } = \e ab Q b (216) 
{M,Q a } = \e ab Q b (217) 
{B,Q a } = a ab Q b . (218) 

5.3 Supermembrane Connection 

The equations for the supersymmetric Chaplygin fluid devolve from a supermembrane La- 
grangian, Lm- We shall give two different derivations of this result, which make use of two 
different parameterizations for the parameterization-invariant membrane action and give rise, 
respectively, to (176) and (|201|). The two derivations follow what has been done in the bosonic 



case in Sections LI and |4 

We work in a light-cone gauge-fixed theory: The supermembrane in 4-dimensional space- 
time is described by coordinates (p = 0,1,2,3), which are decomposed into light-cone 
components = ± X 3 ) and transverse components X 1 {i = 1,2}. These depend on 
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an evolution parameter if = T and two space- like parameters (p r {r = 1,2}. Additionally there 
are two-component, real Grassmann spinors ip, which also depend on r and <p r . In the light- 
cone gauge, X + is identified with r, X~ is renamed 8, and the supermembrane Lagrangian 
is g 

L M = J d 2 <pjC M = - j d 2 <p{VG- ±e rs d r ipad s ip ■ X} (219) 
where G = det G a p; 

g (a„ g„\ = hd T e - (d T xf - 4,d T i> u,\ m 

\G ro g rs J y u r g rs J 

G = gT 

T = 2d T 9 - {d T Xf - 4>d T ?p + g rs uru s 
g r s = d r X ■ d s X , g = detg rs 

u s = d s 9 - lipd s ip - d T X ■ d s X . (221) 

Here d T signifies differentiation with respect to the evolution parameter r, while d r differen- 
tiates with respect to the space-like parameters ip r ; g rs is the inverse of g rs , and the two are 
used to move the (r, s) indices. Note that the dimensionality of the transverse coordinates X 1 
is the same as of the parameters ip r , namely two. 

5.4 Hodographic transformation 



To give our first derivation following the procedure in Section 4.3, we rewrite the Lagrangian 
in canonical, first-order form, with the help of bosonic canonical momenta defined by 

dC M 



dd T X 
dC M 



p = -Ud T X - Uu r d r X (222a) 



n = vVT . (222b) 



dd T 6 

(The Grassmann variables already enter with first-order derivatives.) The supersymmetric 
extension of (|164!) then reads 



C M =p-d T X + ud T e - in^ r v + ^(p 2 + g) + ^e rs d r ^ad s ^ ■ X 

+ u r (p • d r X + Ud r 6 - lU%bd r ^j . (223) 

In (|223j ) u r serves as a Lagrange multiplier enforcing a subsidiary condition on the canonical 
variables, and g = detg rs . The equations that follow from ( |223| ) coincide with the Euler- 
Lagrange equations for ( pl9| ). The theory still possesses an invariance against redefining the 
spatial parameters with a r-dependent function of the parameters. This freedom may be used 
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to set u T to zero and fix IT at —1. Next we introduce the hodographic transformation, as in 



Section 4.3, whereby independent-dependent variables are interchanged, namely we view the 
ip r to be functions of X i . It then follows that the constraint on ( |223| ), which with II = — 1 
reads 



becomes 



p ■ d r X - d r 6 + \i)d r i\) = 



d r X . ( p - V0 + ±vvv> ) = o . 



(224) 



(225) 



Here p, 9 and ip are viewed as functions of X, renamed r, with respect to which acts the 
gradient V. Also we rename p as v, which according to ( |225| ) is 



v = V< 



(226) 



As in Section L2 , from the chain rule and the implicit function theorem it follows that 

8 T = d t + d T X • V (227) 

and according to ( [222a ) (at II = — 1, u r = 0) d T X = p = v. Finally, the measure transforms 
according to d 2 (p — > d 2 r -^=. Thus the Lagrangian for ( ^23[) becomes, after setting u r to zero 
and II to — 1, 

/j2 
— = (v 2 - 9 - v- V9 + \^ + v-V^)- \{v 2 + g) 



\e TS iptf djijj d s x j d r x* 



(228a) 



But e rs d s xW r x i = e 1 ^ detd r x % = e %3 \fg. After y/g is renamed y/2X/p, ( |228a| ) finally reads 



L M = -±=J & 2 r (-p(6 - ±VV0 - \p{Ve - i^V^) 2 - - p - ^^oc x W) . (228b) 

Upon replacing ip by ^(1 — this is seen to reproduce the Lagrange density ( |176| ), apart 
from an overall factor. 



5.5 Light-cone parameterization 

For our second derivation, we return to ( |219D -( |22~1~| ) and use the remaining reparameterization 
freedom to equate the two X 1 variables with the two ip T variables, renaming both as r % . Also 



r is renamed as t. This parallels the method in Section 4.1. Now in ( |2L9| ) ( J22T| ) g ri 
d T X = 0, so that ( |221| ) becomes simply 

G = T = 29-ipip + u 2 (229) 
u = V0 - ±V>VV> • (230) 
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Therefore the supermembrane Lagrangian ( |219| ) reads 

L M = - J d 2 r j^20 -i/nji + (V0 - \^Vi>f + fact x Vt/A . (231) 

Again a replacement of ip by -^=(1 — e)tp demonstrates that the integrand coincides with the 
Lagrange density in (|201| ) (apart from a normalization factor). 



5.6 Further consequences of the supermembrane connection 

Supermembrane dynamics is Poincare invariant in (3+l)-dimensional space-time. This invari- 
ance is hidden by the choice of light-cone parameterization: only the light-cone subgroup of the 
Poincare group is left as a manifest invariance. This is just the (2 + 1) Galileo group generated 
by H, P, M, B, and N. (The light-cone subgroup of the Poincare group is isomorphic to 
the Galileo group in one lower dimension |33].) The Poincare generators not included in the 
above list correspond to Lorentz transformations in the "— " direction. We expect therefore 
that these generators are "dynamical" , that is, hidden and unexpected conserved quantities of 
our supersymmetric Chaplygin gas, similar to the situation with the purely bosonic model. 
One verifies that the following quantities 



D = tH - J d 2 r p6 (232) 

G = j & 2 r{rH - 9Vi - l^aa ■ Vi^) 

= J & 2 r{rH -9T- \ij)cta. ■ Vip) (233) 



are time- independent by virtue of the equations of motion (183)-(186), and they supplement 



the Galileo generators to form the full (3 + 1) Poincare algebra, which becomes the super- 



Poincare algebra once the supersymmetry is taken into account. Evidently (232), (|233j) are 



the supersymmetric generalizations of (105), ( |106| ). 

We see that fluid dynamics can be extended to include Grassmann variables, which also 
enter in a supersymmetry-preserving interaction. Since our construction is based on a super- 
membrane in (3+l)-dimensional space-time, the fluid model is necessarily a planar Chaplygin 
gas. It remains for the future to show how this construction could be generalized to arbitrary 
dimensions and to different interactions. Note that Grassmann Gauss potentials ij) can be 
used even in the absence of supersymmetry. For example, our theory ( |176|) , with the last 
term omitted, posseses a conventional, bosonic Hamiltonian without supersymmetry, while 
the Grassmann variables are hidden in v and occur only in the canonical 1-form. 
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6 One-dimensional Case 



In this section, I shall discuss the nonrelativistic/relativistic models in one spatial dimension. 
Complete integrability has been established for both the Chaplygin gas |53 and the Born-Infeld 



theory [35]. We can now understand this to be a consequence of the complete integrability 
of the Nambu-Goto 1-brane (string) moving on 2-space (plane), which is the antecedent of 
both models. [Therefore, it suffices to discuss only the Chaplygin gas since solutions of the 



Born-Infeld model can then be obtained by the mapping (172)-(173).] 

As remarked previously, in one dimension there is no vorticity, and the nonrelativistic 
velocity v can be presented as a derivative with respect to the single spatial variable of a 
potential 9. Similarly, the relativistic momentum p = v/y/l — v 2 /c 2 is a derivative of a 
potential 6. In both cases the potential is canonically conjugate to the density p governed 
by the canonical 1-form J dx6p. Moreover, it is evident that at the expense of a spatial 
nonlocality, one may replace by its antiderivative, which is p both nonrelativistically and 
relativistically (nonrelativistically p = v), so that in both cases the Lagrangian reads 



L = — ~ J dx dy p(x)e(x — y)p(y) — H . 
For the Chaplygin gas and the Born-Infeld models, H is given respectively by 

^Chaplygin = | dx (l pp 2 + ^ 
^Born-Infeld = j faty ftf + Q 2 ^2 + p 2) _ 

The equations of motion are, respectively 



Chaplygin gas: 



p + -^(pp) = o 

d (p 2 A \ 



d 1 dp 
or — — , + 







Born-Infeld model: 



or 



d_ 

■ 2 dt 



Xe 2 +p 2 




In the above, eqs. ( |239| ) and ( |242[ ) result by determining p in terms of 6 (p = 
and (|24lD, and using that expression for p in (|23?|) and (p40|) . 



(234) 

(235) 
(236) 

(237) 
(238) 
(239) 

(240) 

(241) 
(242) 

from (|238j) 
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6.1 Specific solutions for the Chaplygin gas on a line 

Classes of solutions for a Chaplygin gas in one dimension can be given in closed form. For 
example, to obtain general, time-rescaling-invariant solutions, we make the Ansatz that 9 oc 
1/t. Then (|10S|) or (239) leads to a second-order nonlinear differential equation for the x- 



dependence of 8. Therefore solutions involve two arbitrary constants, one of which fixes the 
origin of x (we suppress it); the other we call k, and take it to be real. The solutions then read 

0(t, x) = --^5- cosh 2 kx . (243) 
2k z t 

[Other solutions can be obtained by relaxing the reality condition on k and/or shifting the 
argument kx by a complex number. In this way one finds that 9 can also be sinh 2 kx, 
sin2 2P7 cos2 but these lead to singular or unphysical forms for p.] The density 
corresponding to (|243|) is found from (|96"|) or Q238]) to be 



The velocity/momentum v = p = jj^d is 



p(t,x) = V2X k }l\ . (244) 
cosh kx 

d 



while the sound speed 



v(t,x) = p(t,x) = — — sinh kx cosh kx (245) 



cosh 2 kx , , . 

s{t,x)= (246) 



is always larger than \v\. Finally, the current j = exhibits a kink profile, 

j(t,x) = -e(t)V2Xtanhkx (247) 
which is suggestive of complete integrability. 



Another particular solution is the Galileo boost of the static profiles (115), ( |116| ): 



p(t, x) = p(x - ut) (248) 



p(t,x) = r^- (249) 
\p — u\ 

Here u is the boosting velocity and p(x — ut) is an arbitrary function of its argument (provided 
p u). Clearly this is a constant profile solution, in linear motion with velocity u. 

Further evidence for complete integrability is found by identifying an infinite number of 
constants of motion. One verifies that the following quantities 



are conserved. 



f / v2A\ n 
lt= / dxp(p± — ) , n = 0,±l,... (250) 
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The combinations p ± are just the velocity (±) the sound speed, and they are known 
as Riemann coordinates. 



6.2 Aside on the integrability of the cubic potential in one 
dimension 

Although it does not belong to the models that we have discussed, the cubic potential for 
1-dimensional motion, V(p) = lp 3 /3, is especially interesting because it is secretly free - a 
fact that is exposed when Riemann coordinates are employed. For this problem these read 
R± = p ± y/2£p and again they are just the velocity (±) the sound speed. In contrast to ( |252j ) 
the Euler and continuity equations for this system decouple: R± = —R±-§^R±- Indeed, it is 
seen that R± satisfy essentially the free Euler equation [compare with ( f42| ) and identify R± 
with v]. Consequently, the solution (f44|)-(|46|) works here as well. 

Recall the previous remark in Section [D] on the Schrodinger group [Galileo © SO(2,l)]: 
in one dimension the cubic potential is invariant against this group of transformations, and in 
all dimensions the free theory is invariant JT^j], |]l8[ . Therefore a natural speculation is that 
the secretly noninteracting nature of the cubic potential in one dimension is a consequence of 
Schrodinger group invar iance. 

Another interesting fact about a one-dimensional nonrelativistic fluid with cubic potential 
is that it also arises in a collective, semiclassical description of nonrelativistic free fermions in 
one dimension, where the cubic potential reproduces fermion repulsion [pq] . In spite of the 
nonlinearity of the fluid model's equations of motion, there is no interaction in the underlying 
fermion dynamics. Thus, the presence of the Schrodinger group and the equivalence to free 
equations for this fluid system is an understandable consequence. 

6.3 General solution for the Chaplygin gas on a line 

The general solution to the Chaplygin gas can be found by linearizing the governing equations 
(continuity and Euler) with the help of a Legendre transform, which also effects a hodographic 
transformation that exchanges the independent variables (t, x) with the dependent ones (p, 9); 
actually instead of p we use the sound speed s = \p2\j p and instead of 9 we use the momentum 




P 

The equations of motion for this system [continuity (|237| ) and Euler ( |238D 1 can be succinctly 
presented in terms of R±: 




Define 








(253) 
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From the Bernoulli equation we know that 



-\p 2 + \s 2 . 



Thus 



tp(p,s) = 6(t,x) + -(p 



2 s 2 ) 



xp 



and the usual Legendre transform rules govern the derivatives. 

dtp 

= tp-x 

op 



dip 
ds 



-ts 



(254) 
(255) 

(256a) 
(256b) 



It remains to incorporate the continuity equation (237) whose content must be recast by 
the hodographic transformation. This is achieved by rewriting equation (237) in terms of 
s = V2X/p: 



ds 



ds 



t s ^p = 

dt dx dx 
Next ( |257| ) is presented as a relation between Jacobians: 

d(s,x) d(t,s) s d(t,p) _ Q 
d(t,x) d(t,x) d(t,x) 

which is true because here dx/dt = dt/dx = 0. Eq. (258a) implies, after multiplication by 

d(t,x)/d(s,p) 

d(t,p) 



(257) 



(258a) 



Q = d(s,x) | d(t,s) _ ^ 

d(s,p) d(s,p) d(s,p) 
dx dt dt 
dp ^ dp ds 



(258b) 



The second equality holds because now we take ds/dp = dp/ds = 0. Finally, from (|255j) , ( |256| ) 
it follows that ( 258b| ) is equivalent to 

d 2 ip d 2 ip 2 dtp 



dp 2 



ds 2 ^ s ds 



. 



(258c) 



This linear equation is solved by two arbitrary functions of p± s (pis being just the Riemann 
coordinates) 



tp(p, s) = F(p + s) - sF'(p + s) + G(p - s) + sG'(p - s) . 



(259) 



In summary to solve the Chaplygin gas equations, we choose two functions F and G, 
construct ip as in (|259| ), and regain s (= y/2X/p), p (= -^-0), and 9 from (|255| ), (|256j ) . In 
particular, the solution (|243|) , (|244|) corresponds to 



F{z) = G(-z) 
where the sign is correlated with the sign of t. 



± — In z 

2k 



(260) 
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53 



6.4 Born-Infeld model on a line 

Since the Born-Infeld system is related to Chaplygin gas by the transformation described in 



Section ^J, there is no need to discuss separately Born-Infeld solutions. Nevertheless, the 
formulation in terms of Riemann coordinates is especially succinct and gives another view on 
the Chaplygin/Born-Infeld relation. 

The Riemann coordinates R± for the Born-Infeld model are contructed by first defining 

1 d „ 

-— — v = p/c = tan (p p 
cox 

a/pc = tainp p (261) 

and 

R± = ip p ± ip p . (262) 



The 1-dimensional version of the equations of motion (|130| ), (|131| ), that is, fl240D , (|241|) can be 
presented as 

d 

R ± = - c ( s inR ) R ± . (263) 
ox 

The relation to the Riemann description of the Chaplygin gas can now be seen in two ways: 
a nonrelativistic limit and an exact transformation. For the former, we note that at large c, 
ip p w p/c, if p a/pc so that 



^Born-Infold ^ \ f _j_ <±\ _ \ ^Chaplygin 
± C\ pJ C ± 



(264) 

A=a 2 /2 



Moreover, the equation ( p63| ) becomes, in view of (264) 



^^Chaplygin _ ^Chaplygin 1 9 ^Chaplygin (265) 

c ± T c dx ± 



so that ( |252| ) is regained. On the other hand, for the exact transformation we define new 
Riemann coordinates in the relativistic, Born-Infeld case by 

K± = csmR± . (266) 

Evidently ( |263j ) implies that 1Z± satisfies the nonrelativistic equations ( |252| ), ( |265| ) when R± 
solves the relativistic equation ( |263j ). Expressing 1Z± and R± in terms of the corresponding 
nonrelativistic and relativistic variables produces a mapping between the two sets. Calling 
PNR; Pnr and pn, /?r the momentum and density of the nonrelativistic and of the relativistic 
theory, respectively, the mapping implied by fl266j ) is 

c 2 prpr 

PNK = 



(p\ + c 2 )(p\c 2 + a 2 ) 
Pnr = i ^(Pi + c 2 )(plc 2 + a 2 ) . (267) 
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As can be checked, this maps the Chaplygin equations into the Born-Infeld equations. But 
the mapping is not canonical. 

We record the infinite number of constants of motion, which put into evidence the (by now 
obvious) complete integrability of the Born-Infeld equations on a line. The following quantities 
are time-independent: 

I± = ac n ~ l [ dx {lfp - ^ , n = 0,±l,... (268) 
J sm (fp cos (fp 



The nonrelativistic limit takes the above into ( 250 ), while expressing J„ in terms of TZ± 



according to ( |266| ) shows that the integrals in (268) are expressible as series in terms of the 
integrals in (|250|). 



In the relativistic model p need not be constrained to be positive (negative p could be 
interpreted as antiparticle density). The transformation p — > —p, p — * —p is a symmetry and 
can be interpreted as charge conjuguation. Further, p and p appear in an equivalent way. As 
a result, this theory enjoys a duality transformation: 

2 

p^±^p p^±—p. (269) 

Under the above, both the canonical structure and the Hamiltonian remain invariant. Solutions 
are mapped in general to new solutions. Note that the nonrelativistic limit is mapped to the 
ultra-relativistic one under the above duality. Self-dual solutions, with p = ±-%P> satisfy 

p = Tc± P (270) 



and are, therefore, the chiral relativistic solutions that were presented at the end of Section 1^2 
In the self-dual case, when p is eliminated from the canonical 1-form and from the Hamiltonian 
with the help of ( |269| ) , one arrives at an action for p, which coincides (apart from irrelevant 
constants) with the self-dual action, constructed some time ago [j37| j 



c 2 



i/d*d*dyp(*M* - v)p(y) - / V?^v^d* 

\ J dtdxdy p{x)e{x - y)p{y) - ^ j dt dx(p 2 (x) + ^) 1 (271) 



6.5 General solution of the Nambu-Goto theory for a (d=l)- 
brane (string) in two spatial dimensions (on a plane) 

The complete integrability of the Chaplygin gas and of the Born-Infeld theory, as well as the 
relationships between the two, derives from the fact that the different models descend by fixing 
in different ways the parameterization invariance of the Nambu-Goto theory for string on a 
plane. At the same time, the equations governing the planar motion of a string can be solved 
completely. Therefore it is instructive to see how the string solution produces this Chaplygin 
solution [Elll. 
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We follow the development in Section 4.3. The Nambu-Goto action reads 



/NG = j 


f dp L N g 




Lng = j 


f dp 1 C NG 




£nG = 


OX^dX^ 

dp* d<pP- 


1/2 



(272a) 
(272b) 
(272c) 



Here X^, fj, = 0, 1,2, are string variables and (pP,p l ) are its parameters. As in Section 4.3 , 
we define light-cone combinations X^ = ± X 2 ), rename X~ as 9, and choose the 

parameterization X + = p° = r. After suppressing the superscripts on p° and X , we construct 
the Nambu-Goto Lagrange density as 

£NG = det v^-^> 2 

u = d^O — drXd^X 
Equations of motion are presented in Hamiltonian form: 



u 



<d v Xf 



V = 



dd-rX 



n = 



NG 



dC 

dd T 6 



(273) 
(274) 

(275) 



d T X = -— p- ud v X 



1 



2n 2 



(p 2 + (d v X) 2 ) - ud„9 



1 



Mf[ d v x 



d T p ■- 
d T U = -d v (uU) 



d^{up) 



and there is the constraint 



pd^x + ad j = o 



(276a) 

(276b) 

(276c) 
(276d) 

(277) 



There still remains the reparameterization freedom of replacing p by an arbitrary function 
of r and p; this freedom may be used to set u = 0, LT = —1. Consequently, in the fully 
parameterized equations of motion Eq. ( 276d| ) disappears; instead of ( 276a ) and ( 276c ), we 
have d T X = p, d T p = d 2 X, which imply 



(d 2 - d 2 )X = 



(278a) 



( |276b|) reduces to 



d T e = \ [{d T xf + {d^xf 



(278b) 
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and the constraint (j277|) requires 



d v = d T Xd v X . 



(278c) 



(279) 



(280) 



This completes the description of a string moving on a plane. But we need to convert this 



Solution to (278a) is immediate in terms of two functions F±, 

x(t, cp) = F + (t + ip) + F_(r - <p) 

and then ( g7gg) , (^78^ ) fix 6: 

9(r,<p) = / dz[F' + (z)] 2 + / d*[Fi( 



information into a solution of the Chaplygin gas, and we know from Section 4.3 that this can 
be accomplished by a hodographic transformation: instead of X and 9 as a function of r and 
(p, we seek <p as a function of r and X, and this renders 6 to be a function of r and X as well. 
The density p is determined by the Jacobian \dX/d<p\. 

Replace r by t and X by x and define </? to be f(t, x). Then from ( p79[ ) it follows that 



(281) 



x = F+(t + /(t,x)) +F_(t-/(t,s)) . 

This equation may be differentiated with respect to t and x, whereupon one finds 

F'(t + f) + FL(t-f) 



df 
dt 

df 

dx F' + {t + f)-F'_{t-f) 



F' + {t + f)-F'_{t-f) 
1 



(282a) 
(282b) 



Thus the procedure for constructing a Chaplygin gas solution is to choose two functions 
F±, solve the differential equations (|282j) for /, and then the fluid variables are 



6(t,x) 
V2X 



P 



t+f(t,x) rt-f(t,x) 

[*+(*)] 2 d*+ / [F'-{z)] 2 dz 



\F'(t + <p)-FL(t-<p)\ 



(283) 
(284) 



One may verify directly that ( p83| ) and ( |284| ) satisfy the required equations: Upon differ- 
entiating (|283j ) with respect to t and x, we find 



89 
di 

09 



df 



-2FiFi 



(285a) 
(285b) 
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The second equalities follow with the help of fl282|). From (|285|) one sees that 



89 1 {d6\2 



1 1 I?' 



X 



dx) 



(286) 



the last equality being the definition ( |284| ) . Thus the Bernoulli (Euler) equation holds. For 
the continuity equation, we first find from ( 284 ) and (|285| ) 

dp 



d 


V2X 


dtF\ 


-FL 


V 


T\ 


(K- 


-F'_) 2 


2\ 


/2A 



Fl 1 + 



df 

at 

F'iF'_ + F'LFL 



Fill 



at) 



(287a) 



d t d8\ d 
dxVdx) dx\ 

= =F 



(±V2X 



F' + +F'_ 
F'-F' 



- A ( p// p/ I p/ \ ®f 
(F'-FL) 2[+ - + - +) dx 



(287b) 



The last equalities follow from ( p82| ) ; since ( 287a| ) and ( 287b| ) sum to zero, the continuity 
equation holds. 

We observe that the differentiated functions F± are just the Riemann coordinates: from 
( 285b ) and (|284| ) [with the absolute value ignored] we have 



p ± 



P 



R± = 2FL 



± ■ 



(288) 



Also it is seen with the help of ( |282| ) that the Riemann formulation ( 252 ) of the Chaplygin 
equations is satisfied by 2F'±. 

The constants of motion (|250|) become proportional to 



It oc / dx 



1 



F' - F' 



~x 



Finally we remark that the solution (|243| ) , ( p44| ) corresponds to 

„ , n -w—t / \ In z 
F+(z) = -F-(z)=±— . 



(289) 



(290) 
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There exists a relation between the two functions F and G in (|259|), which encode the Chaplygin 



gas solution in the linearization approach of Section 6.3, and the above two functions F±, which 
do the same job in the Nambu-Goto approach. The relation is that 2F' + is inverse to 2F" and 
2F'_ is inverse to 2G" , that is, 



2F"[2F' + {z)]= z 
2G"[2F'_{z)]= z 



(291) 



Problem 7 Derive ( p9l|) . Verify this relation with (|260|) and (|29Cl ). 
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7 Towards a Non-Abelian Fluid Mechanics 

Fluid mechanics and fluid magnetohydrodynamics may very well describe the long-wavelength 
degrees of freedom in a quark-gluon plasma. Moreover it is plausible that the group (color) 
degrees of freedom remain distinct in that regime, so that one should incorporate them in 
the fluid approximation. In this way one is led to think about constructing non-Abelian fluid 
mechanics and (color) magnetohydrodynamics. In this section we describe an approach to this 



task [38]. In the course of development, we encounter and solve an interesting mathematical 
problem: how to parameterize a non-Abelian vector potential so that the non-Abelian Chern- 
Simons density becomes a total derivative, and the volume-integrated Chern-Simons term is 
given by a surface integral. Obviously this is the non-Abelian generalization of the similar 
Abelian problem, which is solved by presenting the Abelian vector potential in Clebsch form. 
So we shall determine the non-Abelian version of the Clebsch parameterization. 



7.1 Proposal for non-Abelian fluid mechanics 

We review our Lagrange density for relativistic Abelian fluid mechanics (|144|): 



C = -j% - f(^F%) . (292) 

The equation of state is encoded in the function /. For free fluid motion f(y/J^J^) = c^j^j^. 
Here is the matter current and a M is an auxilliary 4- vector, which is presented in the form 

a M = d^6 + adpP ■ (293) 



The time component do, involving time derivatives, determines the canonical 1-form; the 
spatial components a are in the Clebsch parameterization, as is needed for overcoming the 
obstacle created by a Casimir invariant of the fluid in the algebra (|14|), ([l6|). Another way of 
characterizing the parameterization of the vector a is that it casts the Chern-Simons density 
of a, namely, a • V x a into total derivative form: V# • (Va x V/3) = V • (0V x a). 

For a non-Abelian generalization, it is plausible to suppose that the current 4-vector ac- 
quires an internal symmetry index: J a ; correspondingly, the auxiliary 4-vector must also 
acquire an internal symmetry index: A". It remains to give a rule for parameterizing A®, 
which generalizes the Abelian rule ( p93| ) . 

Our proposal - and it is a speculative one, since at this stage we have no derivation 
from microscopic considerations - is that should be written in a form so that its non- 
Abelian Chern-Simons density, CS(A) = A a dA a +\f abc A a A h A c , is a total derivative (f abc are 
the structure constants of the group). This leads us to the purely mathematical problem of 
constructing a parameterization for a non-Abelian vector that ensures this property. 
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7.2 Non-Abelian Clebsch parameterization 

(or, casting the non-Abelian Chern-Simons density into total derivative form) 

We enquire whether it is possible to parameterize the non-Abelian 1-form, A a , such that the 
Chern-Simons 3- form is a total derivative (is exact): 

CS(,4) = A a dA a +\f abc A a A h A c = dtt . (294) 

That this should be possible follows from the observation that the left side of (|294|) is a 3-form 
on 3-space; hence it is closed, because a 4-form does not exist in 3-space. [Of course on a 4- 



dimensional space the exterior derivative of ( 294 ) is proportional to the non-Abelian anomaly 



(Chern-Pontryagin density) But a closed form is also exact, at least locally; this justifies 

the right side of (|294|) , 



How this works in the Abelian case has already been explored in Sections 2.4 and 2.5: the 
Clebsch parameterization (|58|), ( |60| ) for A leads to the desired result. But the generalization 
of (|58|), (|60|) for a non-Abelian 1-form is not evident. However, at the end of Section 2J3, an 



alternative approach is presented, wherein the Abelian 1-form is projected from a non-Abelian 
pure gauge 1-form. This construction can be generalized to the non-Abelian case and yields 
the sought-for parameterization. 

The mathematical problem can therefore be formulated in the following way: For a given 
group H, how can one construct a potential A^ = (Aq,A?) such that the non-Abelian Chern- 
Simons integrand CS(A) is a total derivative? Here we shall only sketch the solution to the 
problem, referring those interested to Ref. |^] for a detailed discussion. 

In the solution that we present, the "total derivative" form for the Chern-Simons density of 
A a is achieved in two steps. The parameterization, which we find, directly leads to an Abelian 
form of the Chern-Simons density: 

A a dA a + \f abc A a A b A c = 7 d 7 (295) 



for some 7. Then Darboux's theorem |10| (or usual fluid dynamical theory |Tl[]) ensures that 
7 can be presented in Clebsch form, so that 7 d7 is explicitly a total derivative. 

We begin with a pure gauge g~ l dg in some non-Abelian group G (called the Ur-group) 
whose Chern-Simons integral coincides with the winding number of g. 

W(g) = J d 3 r CS^ 1 dg) = ^J trig" 1 dgf (296) 

We consider a normal subgroup H C G, with generators I a , and construct a non-Abelian 
gauge potential for H by projection: 

A a oc tr(/ a 9 ~ 1 dg) . (297) 

Within H, this is not a pure gauge. We determine the group structure that ensures the 
Chern-Simons 3-form of A a to be proportional to tr(g -1 dg) 3 . Consequently, the constructed 
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non-Abelian gauge fields, belonging to the group H, carry quantized Chern-Simons number. 
Moreover, we describe the properties of the Ur-group G that guarantee that the projected 
potential A a enjoys sufficient generality to represent an arbitrary potential in H. 

Since tr^ -1 dg) 3 is a total derivative for an arbitrary group (although this fact cannot in 
general be expressed in finite terms |39f|) our construction ensures that the form of A a , which is 



achieved through the projection ( 297 ), produces a "total derivative" expression (in the limited 
sense indicated above) for its Chern-Simons density. 

Conditions on the Ur-group G, which we take to be compact and semi-simple, are the 
following. First of all G has to be so chosen that it has sufficient number of parameters to 
make ti{I a g~ l dg) a generic potential for H. Since we are in three dimensions, an //-potential 
Af has 3 x dim// independent functions; so a minimal requirement will be 

dim G > 3 dim H . (298) 

Secondly we require that the //-Chern-Simons form for A a should coincide with that of g~ l dg. 
As we shall show in a moment, this is achieved if G/H is a symmetric space. In this case, if 
we split the Lie algebra of G into the //-subalgebra spanned by I a , a = 1, . . . , dim H, and the 
orthogonal complement spanned by S A , A = 1, . . . , (dim G — dim H), the commutation rules 
are of the form 

[I* , S A ] = h aAB S B (299b) 
[S A , S B ] = N h aAB I a . (299c) 

(h a ) AB form a (possibly reducible) representation of the //-generators I a . The constant N 
depends on normalizations. More explicitly, if the structure constants for the Ur-group G 
are named f abc , a,b,c = 1, . . . , dimG, then the conditions ( |299a -c) require that f abc vanishes 
whenever an odd number of indices belongs to the orthogonal complement labeled by A, B, ... 
Moreover, f abc are taken to be the conventional structure constants for H and this may render 
them proportional to (rather than equal to) f abc . 
We define the traces of the generators by 

tr(/ a / 6 ) = -N^ , tr(S A S B ) = -N 2 5 AB 

tr(I a S A ) = . (300) 

We can evaluate the quantity tr[S A , S B ]I a = tr S A [S B , I a ) using the commutation rules. This 
immediately gives the relation N\N = N 2 . 

Expanding g~ 1 dg in terms of generators, we write 



9 



- 1 dg = (I a A a + S A a A ) (301) 



which defines the //-potential A a . Equivalently 



A a = -^-ti(I a g- l dg) (302) 
iVi 
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From d(g 1 dg) = —g 1 dg g 1 dg, we get the Maurer-Cartan relations 

F a = dA a +y abc A b A c = -—h aAB a A a B 

da A +h aBA A a a B = . (303) 

Using these results, the following chain of equations shows that the Chern-Simons 3-form for 
the i/-gauge group is proportional to tr(<7 _1 dg) 3 : 

1 : {A a dA a +\f abc A a A b A c ) = — l -^(A a dA a +2 A a F a ) 



16tt 2 " iJ ' 48?r 



' (A a dA a -Nh aAB A a a A a B ) 



487T 2 



48?r 



' (A a dA a +N da A a A ) 



1 N 
48vrn^ tr(j4dj4) + ^ tr(daa) 

--^-tr(AdA+ada) 
.-J^tr g-ldg d^ 1 dg) 

' trig' 1 ■ (304) 



48vr 2 iVi 



In the above sequence of manipulations, we have used the Maurer-Cartan relations ( |303j ) , 
which rely on the symmetric space structure of ( 299a| -c), and the trace relations ( |300| ) , along 
with NiN = N 2 . 

We thus see that J CS(A) is indeed the winding number of the configuration g G G. Since 
tr(g -1 dg) 3 is a total derivative locally on G, the potential ( |302| ), with the symmetric space 
structure of ( 299a| -c), does indeed fulfill the requirement of making CS(j4) a total derivative. 



It is therefore appropriate to call our construction ( |302D a "non-Abelian Clebsch parameteri- 
zation" . 

In explicit realizations, given a gauge group of interest H, we need to choose a group G 
such that the conditions ( |298| ), ( |299a| -c) hold. In general this is not possible. However, one can 
proceed recursively. Let us suppose that the desired result has been established for a group, 
which we call H2- Then we form H C G obeying (299a-c) as H = H\ x H2, where Hi is the 



gauge group of interest, satisfying dimG > 3 dimH\. For this choice of H, the result ( |30 
becomes 

CS(^) + CS^) = ^-L^- trGT 1 dg) 3 (305) 



But since CS(i?2) is already known to be a total derivative, ( 305 ) shows the desired result: 
CS(-Hi) is a total derivative. 

To see explicitly how this works we work out the representation for a SU(2) ~ 0(3) potential 
Af, which possesses nine independent functions. 
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We take G = 0(5), H = 0(3) x 0(2). We consider the 4-dimensional spinorial representa- 
tion of 0(5). With the generators normalized by tr(t a t b ) = — S ab , the Lie algebra generators 
of 0(5) are given by 



I a = 





f a a 







\ 


a a 













1 


1 


( o 





(306) 

S A = — — ( u u ] S A = — !— ( aA | 

iV2\a A J iV2 \ J 

cr's are the 2x2 Pauli matrices. I a generate 0(3), with the conventional structure constants 
e abc , and 1° is the generator of 0(2). S, S are the coset generators. 

A general group element in 0(5) can be written in the form g = Mhk where h £ 0(3), 
k £ 0(2), and 

-, / 1 — |(w x w) ■ a —w-a \ 

M= (307) 
J I + w • w - |(w x w) 2 y w ■ a 1 + |(w x w) • a J 

w a is a complex 3-dimensional vector, with the bar denoting complex conjugation, w • w = 
w a w a and (w x w) a = e abc w b w c . The general 0(5) group element contains ten independent 
real functions. These are collected as six from M (in the three complex functions w a ), three 
in h, and one in k. 

The 0(3) gauge potential given by — tr(/ a <7 _1 dg) reads 

A a = R ab (h) a b + (h- 1 dh) a 

a 1 j w a dw •(wxw)+Mi° dw -(w x w) 

1 + w • w - |(w x w) 2 ) 2 

+ e abc {dw b w c -w b dw c ) \ 



(308) 



where R ab (h) is defined by hl a h~ x = R ab h b and k does not contribute. A a is the /i-gauge 
transform of a a , which depends on six real parameters (w a ). The three gauge parameters of 
h 6 0(3), along with the six, give the nine functions needed to parameterize a general 0(3)- 
[or SU(2)-] potential in three dimensions. The Chern-Simons form is 



CS(^) = -{A a dA a +\e abc A a A b A c ) 

L : {a a da a +\e abc a a a b a c )-d 



' (dhh- 1 )*^) 



16vr 2 



16vr 2 v ' z L 

+ ^tr(h- 1 dhf (309) 
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The second equality reflects the usual response of the Chern-Simons density to gauge trans- 
formations. Using the explicit form of a a as given in ( |30§| ), we can further reduce this. Indeed 
we find that 

(w x dw) • p + (w x dw) • p 



a a da a +\e abc a a a h a c 



(-2 



Pk 



Defining an Abelian potential 



[1 + w • w 4 

\eijk dw* dw 7 

w • dw — w • dw 

1 + w • w — l(w x w) 2 



i^W X w) 2 l 2 



we can easily check that a da reproduces ( plCf ). In other words 



CS(A) 



1 



16vr 2 



a da + d 



(dhh 



16vr 2 



+ 



1 



48vr 2 



tr^d/i) 2 



(310) 



(311) 



(312) 



If desired, the Abelian potential a can now be written in the Clebsch form making a da into 
a total derivative, while the remaining two terms already are total derivatives, though in a 
"hidden" form for the last expression. This completes our construction. 



7.3 Proposal for non-Abelian magnetohydrodynamics 

We return to our construction of a Lagrange density for non-Abelian kinetic theory. As 
explained in Section 7.1, a plausible non-Abelian generalization for fl292|) is 

£ = _ r l tr /Y^- c yji^ (313) 
iVi 

where for simplicity we have taken the "free" form for /. When the desired group is SU(2), g 
is an 0(5) group element, as detailed in Section 7.2. 

Magnetohydrodynamics is achieved by introducing a further interaction with a dynamical 
gauge potential A^. This is accomplished by promoting the derivative of g to a gauge-covariant 
derivative, gauged on the right 



1 



magnetohydrodynamics — tr(/ g D^g^ C-y J^ a J^ ^ J~ uu (314) 

with 

D^g = d^g + egAp ■ (315) 

A^ = A a ^I a are independent, dynamical gauge potentials (not given by g) leading to the field 
strengths T^ v . The gauge transformation properties by the gauge function h are 

g > = g h A' = h~ l Ah + -h- 1 dh 

e 

J'*P = h^J^Ph . (316) 
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We expect that the Lagrangian ( |314| ) will describe non-Abelian magnetohydrodynamics, 
namely the dynamics of a fluid with non-Abelian charge coupled to non-Abelian fields. [The 
Abelian version of (|314| ) does indeed describe ordinary magnetohydrodynamics.] This gluon 
hydrodynamics can be useful for non-Abelian plasmas such as the quark-gluon plasma. Details 
of ( |314| ) and possible applications are under further study. 
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Solutions to Problems 

Problem | The imaginary part of the Schrodinger equation gives the continuity equa- 
tion in the form p + V • (pV#) = 0. This identifies the velocity v as V0, that is, v is 
irrotational and there is no vorticity. The real part becomes the Bernoulli equation 9 + 
7j(V#) 2 = ^-p _1 / 2 V 2 p 1 / 2 , whose gradient gives the Euler equation and identifies the force / 

as V( ^ /9 -l/2 V 2 /9 l/2 ) 

Problem § j = pv with v = V0. 

Problem |3| £schr6din g er = Op - \p{V9) 2 - ^ Vp-Vp w ] iere trie time derivative of i\ - p9 
has been dropped. 

The results in the solutions to Problems [l] and |3| are called the Madelung formulation of 
quantum mechanics |4C]] . 

Problem | CS(A) = e ijk di^dj cos @d k h{r) 

(a) Extracting the first derivative leaves CS(A) = diV*, V* = e l ^ k ^dj cos Qd k h(r). (This is 
true because e^ k didj cos G = = e l ^ k did k h(r), since cos0 and h(r) are nonsingular.) 
Note that V* = £ ier $(-^)ti (r) = <P*$(± sin 9)/i(r). Since V r a = 0, the surface 
integral does not contribute. However, since $ is multivalued, there is a contribution from 
the $ integral: / d 3 r CS(A) = r 2 dr /* sin 9 d6 J Q 2w d^^^^S) x (i sin 8)/i'(r) = 
4ir[h{R) - h{0)]. 

(b) Extracting the second derivative leaves CS(A) = djVjj — e lJ (c^O) cos Qd k h(r). The 
last term is present, owing to the singularity of at the origin, which gives rise to 
e ki Widj<f> = 5 k3 2TT5(x)6{y). (See |IJ.) Also we have VJ* = ^' fc di$cosed fc /i(r) = 
e* jr ( — ETe) cosB/i'(r) = — ^5^ e cot Qh'(r). Again there is no r-component to contribute 
to the surface integral, but the second, singular term leaves 

J d 3 rCS(A) = J dV(27r)5(x)<5(y) cos JU(r) 
f R z 8 

= 4vr f dz ^-h{z) = 4vr[/i( J R) - h{0)] . 
Jo oz 

(c) Extracting the last derivative leaves 

CS(A) = 8 k V k - e ijk (d k di<S>)dj cos 6%) , 
V k = £ ijk d .^ d . cog Qh ^ 

= e* dk ( -J- )(-isine)%) 
V r sin H / \ r / 
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Here the surface integral contributes Airh{R). The singular term is 

rs 

-5 j3 (27T)5(x)5(y)(d j cos @)h(r) = -(2tt)<5(z)%) (— ^(M) 

= -(4ir)5 3 (r)h(0) . 

Hence this contribution to the spatial integral is — 4nh(0), for a total of &ir\h(R) — /i(0)] . 

Problem [5] In the Clebsch parameterization, B = Va x V/3, and 5 A = VJ7 + + 
Therefore 

B 5 A = B V5 7 + B ■ aV5(3 + (B • V(3)5a 

= V • (B5j) + V • (Ba60) - (B • Va)5(3 + (B • V ' (5)5a . 

The last two terms vanish, so J d 3 r A ■ B is the surface term J dS -B(5j + a.80) with no 
contribution from the bulk (finite r space). This of course is consistent with the Chern- 
Simons integral being a surface term, since 6^ f d 3 r A ■ B = J d 3 r B ■ 5 A. When demanding 
the variation of | J d 3 r B 2 to vanish, we first find J* d 3 r(VxB)-i5A = 0. When 5 A is arbitrary, 
this condition implies the vanishing of V x B. However, in the Clebsch parameterization, all 
we can conclude is that VxBis proportional to B, with a position-dependent proportionality 
factor VxB = fiB. Taking the divergence shows that B ■ V/i = 0, that is, n can be a function 



of the magnetic surfaces; see (p4[), (69). 



Problem § 6(t,r) = r 2 /2t. Time rescaling: w {t,r) = e w 0(T,r), T = e w t; e w 9(T,r) = 
e?r 2 /2ePt = r 2 /2t = 6(t,r). 

Space-time mixing: u (t, r) = 0(T, R),T = t + u r + \uj 2 6(T, R)=t + u-r + uj 2 R 2 /AT, 
R = r + u>9(T, R) = r + u>R 2 /2T. Squaring the second equation gives R 2 = r 2 + u> ■ rR 2 /T + 
lu 2 R 4 /AT 2 . Multiplying the first equation by R 2 /T gives R 2 = tR 2 /T + u -rR 2 /T + oj 2 R i /AT 2 . 
Comparing the two shows that R 2 /T = r 2 /t or 6(T, R) = 0(t, r). 



Problem ^ From ( p56|) and (|259| ) we learn that t = F" + G" , x = (p + s)F" - F' + 



(p — s)G" — G' , where F is a function of p + s = 12+ and G is a function of p — s = 12_. 
Differentiating these equations with respect to t and x, it follows that 1 = F"'R + + G"'R-, 
= F"'®R± + G'"^, = R + F"'R + + R_G"'R- , 1 = R+F'"^- + R_G'"^, which in turn 
imply = 1/(12+ - 12_)F"', ^ = -1/(12+ - 12_)G"', and R± = -R T ^ [the Riemann 
equation (|252|) again]. 

On the other hand, the functions F±(t db /) describing the Chaplygin gas solution from 
the Nambu-Goto equation are related to R± by ( p88| ): 12+ = 2F±. Hence ^ = ±2F^|| = 
±2F'l/{F' + -FL) = ±4F£/(R + -R_), where ( gggg ) is used. It follows that 4F + »F'"(2F + (z)) = 
i(2F"(2F|(^))) = 1 and 41^(z)G'"(21^)) = £ (2G"(2Fi(z))) = 1, or 2F" (*+(*)) = 2 
and 2G"(2Fi(z)) = z. When F+(z) = lnz/2A;, 2F' + {z) = 1/zk; with F(z) = Jrhiz, 
2F"(z) = and 2F" (2F' + {z)) = z; similarly for F-(z) and G{z). 
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